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We introduce the pseudohole and heavy-pseudoparticle operator algebra that generates all 
Hubbard-chain eigenstates from a single reference vacuum. In addition to the pseudoholes already 
introduced for the description of the low-energy physics, this involves the heavy pseudoparticles asso- 
ciated with Hamiltonian eigenstates whose energy spectrum has a gap relatively to the many-electron 
ground state. We introduce a generalized pseudoparticle perturbation theory which describes the 
relevant finite-energy ground state transitions. In the present basis these excitations refer to a small 
density of excited pseudoparticles. Our operator basis goes beyond the Bethe-ansatz solution and it 
is the suitable and correct starting point for the study of the finite-frequency properties, which are of 
great relevance for the understanding of the unusual spectral properties detected in low-dimensional 
novel materials. 

PACS numbers: 72.15. Nj, 72.90.+y, 05.30. Fk, 03.65. Ca 



I. INTRODUCTION 



O , 

The unusual spectral properties of one-dimensional electronic quantum liquids imply that they cannot be 
described by one-electron, Fermi-liquid-like models Further, descriptions of these quantum systems in terms 

£NJ \ of exotic excitations such as holons and spinons , or pseudoparticles |?],|| |] , have been limited to the low-energy 

■ Hilbert subspace, and both bosonization |l],|| and conformal-field theory techniques ]To| , pl] , ^2|JT^ ] also apply only in 
that limit. However, low-energy studies cannot describe the finite-frequency properties, which are of great interest 
for the understanding of the unusual properties detected in real low-dimensional novel materials. Thus developing a 
microscopic operator description of these systems for arbitrary energies remains an important theoretical challenge. 

For more than sixty years the "Bethe ansatz" (B A) |l4| has played a central role in the analytic solution of a variety of 
"integrable" many-body problems in condensed matter physics and quantum field theory [ p^l|Jl^Jl^ , [r^ , |l8|] . Techniques 
Q\ , based on the exact BA solutions were used to describe excited states, thermodynamics, and low-energy correlation 

■ functions [p^ , p^|Jr^Jl9[| and remain an active subject of study today. However, the BA solution only provides limited 
information on the many-particle problem. For instance, it does neither provide off-diagonal operator expressions nor 

H ■ matrix elements. 

In this paper we introduce a suitable pseudohole and heavy-pseudoparticle operator algebra which generates all 



Hubbard-chain |18 1£|] eigenstates from a single pseudohole and heavy-pseudoparticle vacuum. This is the pseudohole 
vacuum already introduced in Ref. p0[ | for the low-energy Hilbert subspace spanned by the lowest-weight states 
(LWS's) and highest-weight states (HWS's) of the SU{2) ?7-spin and spin algebras pl| , p2| , p3 24 of type I. (See below 



the definition of types I and II LWS's and HWS's.) That vacuum is the half- filling and zero-magnetic-field ground 
state (GS). 

Recent investigations have established that although not themselves complete, the Hamiltonian eigenstates which 
refer to the B A solution can be extended to form a complete set of states for the one-dimensional Hubbard model |23| . 
However, the vital and to date open issue of how to describe specifically those excitations that dominate response 
and transport at finite energies in integrable quantum liquids and related models remains an open problem. Our 
generalized operator algebra is complete and allows the construction of a pseudoparticle perturbation theory which 
is the correct starting point to solve this problem. The main point is that the GS transitions which determine and 
dominate response and transport at finite energies involve a small density of excited pseudoparticles only. Therefore, 
the description of these relevant transitions is within the range of our perturbation theory. In this paper we introduce 
the generalized operator representation. The application of our approach to study the finite-energy spectral and 
transport properties of the Hubbard chain will be presented elsewhere [^5|,^6| . 

In contrast to three-dimensional Fermi liquids the low-energy excitations of one-dimensional many-electron 

quantum problems are at zero magnetic field characterized by charge - spin separation ^,^,^|J6| . This can be interpre- 
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tated in terms of holon and spinon modes flplp^] . The BA solution of the Hubbard chain |T^] at zero magnetic field 
and chemical potential has allowed the identification and study of the holon and spinon excitations and corresponding 
symmetry transformations ||,|2(J. (The spinon excitations of the Hubbard chain are similar to the corresponding 
spinon excitations of the spin 1/2-isotropic Heisenberg chain p7f.) 

On the other hand, at finite magnetic fields the charge and spin separation is replaced by a more exotic c and s 
separation ||. Here c and s refer to orthogonal small-momentum and low-energy modes which couple both to the 
charge and spin channels. These can be described in terms of pseudoparticle and pseudohole excitations (2(]]. In Ref. 
pof it was shown that the zero-magnetic field holon and spinon excitations are limiting cases of the general pseudohole 
excitations. This has allowed to obtain the exact expression of the electron operator in terms of holon and spinon 
operators pp| . This study has also provided the relation between Hamiltonian symmetry and the transformation of 
the elementary excitations at finite magnetic fields. (At finite magnetic fields the usual holon and spinon picture 
breaks down j2(j].) 

The present pseudohole and heavy-pseudoparticle operator algebra is a generalization of the low-energy pseudopar- 
ticle operator representation of Refs. |7|J||| and pseudohole representation of Ref. [^0). For the general case of the 
Hubbard chain at finite magnetic field and chemical potential both the non-LWS's and non-HWS's of the ?7-spin and 
spin algebras have energy gaps relative to the corresponding canonical-ensemble GS p0| , |2"8| ] . (The expressions of these 
gaps are determined in this paper.) The LWS's and HWS's of these algebras can be classified into two types, the 
states I (or LWS's I and HWS's I) and the states II (or LWS's II and HWS's II). While the Hamiltonian eigenstates 
I are described only by real BA rapidities, all or some of the rapidities associated with the eigenstates II are complex 
and non real. Since for finite magnetic field and chemical potential the states II have energy gaps relative to the 
corresponding canonical-ensemble GS, in that case the low-energy physics is exclusively determined by the states I 
p0| . Our algebra refers to the whole parameter space and is expressed most naturally in terms of operators describ- 
ing the pseudohole excitations introduced in Ref. |2C| ], the new heavy-pseudoparticle excitations, and the associate 
topological momentum-shift excitations. 

We demonstrate that not only the familiar BA solution of the Hubbard chain |lj,|l^,^] can be interpreted naturally 
at all energy scales in terms of the pseudohole and heavy-psepseudoparticle basis, but that such basis goes beyond that 
solution. We then establish that the generalized pseudoparticle perturbation theory, which is also naturally described 
in terms of the new operator algebra, can be used to study the low (w — luq) energy excitations. Here ujq represents 
the energy gaps of particular finite-energy states II and non-LWS's and non-HWS's which proven to dominate the 
finite-energy response functions f25^,|26| . 

In contrast to the states I which at finite values of the chemical potential and magnetic field span the low-energy 
Hilbert subspace ^0|, the LWS's and HWS's II are characterized by finite occupancies of the heavy-pseudoparticle 
bands, whose energy spectrum has a finite gap. In addition, our pseudohole and heavy-pseudoparticle operator basis 
generates the non-LWS's and non-HWS's of the r;-spin and spin algebras absent in the BA solution. It goes beyond the 
BA solution jl8|,|l9|,^3j because it i)- generates from a single vacuum all Hamiltonian eigenstates (including the above 
ry-spin and spin non-LWS's and non-HWS's); ii)- provides relevant information on non-diagonal operators [g^,|6); and 
iii)- since the relevant finite-energy transitions involve a small density of pseudoparticles, it allows the direct study of 
the finite-energy transport properties p5| , |26| j . 

Importantly, we find that at all energy scales there is only zero-momentum pseudoparticle forward scattering in the 
Hubbard chain. Although this exact property holds only for integrable systems, we believe that our pseudohole and 
heavy-pseudoparticle operator description will provide the most suitable starting point for description of the unusual 
finite-frequency properties of both integrable and non-integrable electronic Luttinger liquids. In the case of the latter 
liquids this could be achieved by means of the bosonization of the pseudoparticle theory Q . 

We fully define the pseudohole and heavy-pseudoparticle statistics by evaluating their infinite-dimensional 
statistical- interaction matrix p9| . They are neither fermions nor bosons, yet they obey an anticommuting alge- 
bra. Furthermore, the introduction of topological-momentum-shift operators |2(],|3(J allows the use of a suitable 
second-quantization formalism. 

Since the Hamiltonian-eigenstate transitions involve only our pseudohole, heavy-pseudoparticle, and topological- 
momentum-shift operators, all operators acting onto the Hilbert space can be written in terms of the former operators. 
However, the BA alone does not provide directly the expressions of off-diagonal [in the energy eigenstate basis] part of 
the physical operators. Among the specific results we quote in this paper are the expressions in terms of pseudohole 
and heavy-pseudoparticle operators for the Hamiltonian, momentum operator, and six generators of the ?7-spin and 
spin 5*0(4) algebra (yet four of these generators are off diagonal in the energy representation). Further, the algebra 
can provide diagonal part of all operators acting in the Hilbert space. On the other hand, combining BA with the 
pseudoparticle operator basis and the invariances of the quantum liquid around particular values of momentum and 
energy provides important information on correlation functions p5| , p6[ . (This goes beyond BA because these functions 
are determined by the off-diagonal part of the physical operators.) Our pseudohole and heavy-pseudoparticle operator 
representation is the natural generalization of the low-energy pseudoparticle description which was studied in detail 
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in Refs. 0,||@||,|5|j36 1 . 

In Section II we introduce the pseudohole and heavy-pseudoparticle operator description. In this operator represen- 
tation all Hamiltonian eigenstates have a simple Slater-determinant like form of pseudohole and heavy-pseudoparticle 
levels. 

In Section III we introduce the key concept of generalized ground state (GGS) of which the GS is a particular case. 
We also define the pseudohole and heavy-pseudoparticle statistics J2^] . 

In Section IV we write the Hamiltonian in the pseudohole and heavy-pseudoparticle operator basis. For this we 
relate the rapidity numbers used in the BA |]l9| to the pseudomomentum-distribution operators and use this relation 
to write the Hamiltonian in terms of the pseudoparticle operators. We generalize the rapidity operators for states out 
of the BA solution and show that the same rapidity eigenvalue describes a SU(2) r;-spin tower or a SU(2) spin tower 
of Hamiltonian eigenstates. 

In Sec. V we introduce the generalized pseudoparticle perturbation theory for Hilbcrt subspaces of arbitrary energy. 
This includes evaluation of the heavy-pseudoparticle energy gaps. We discuss the effects of normal-ordering (relative 
to the suitable GGS or GS) and show that, in terms of the pseudoparticle operators, the normal-ordered Hamiltonian 
has in each subspace a simple form. This permits the introduction of a systematic, non-singular pseudoparticle 
perturbation theory valid at low (w — u>o) energy and establishes directly the universal form of the Hamiltonian for 
integrable electronic quantum liquids. 

Finally, in Sec. VI we present the discussion and concluding remarks. 



II. THE GENERALIZED PSEUDOHOLE AND HEAVY-PSEUDOHOLE OPERATOR BASIS 

We consider the Hubbard chain Jl8|,|l9| in a magnetic field H and chemical potential /i 

H = H SO (4) + 2^S C Z + 2^HS S Z , (1) 

where 

1„ + 1, 



HsO{4) 



3,v 



In the absence of the chemical-potential and magnetic-field terms the Hamiltonian (1) reduces to (2) and has 
50(4) = SU{2) <g) SU(2)/Z 2 symmetry f|||f|i|j|§ . In equation (1) we used the notations rj = S c and 5 = S s for 
ry-spin and spin, respectively. With this notation the generators of the 77-spin and spin 50(4) algebra are |2l|p^ , p3| ] , 

SI = ~\[N a S c _= E< l^'y:, , S% = EC" 1 ) 3 '^^ ' ( 3 ) 

" 3 3 

for 77-spin and 



SI 



~5>#„ S!=£ c Jt<*i> ^ = E c i^T, (4) 



for spin. Since N a is even, the operator S c z + S s z has only integer eigenvalues and all half-odd integer representations 
of SU(2) (gi 5/7(2) are projected out |2l"| , |22| , p3| . This selection rule excludes transitions with AS Z + AS Z half-odd 
integer. 

In equations (2) — (4) the operator c - a and Cj a creates and annihilates, respectively, one electron of spin projection 
a (a refers to the spin projections =\ , J, when used as an operator or function index and is given by a = ±1 
otherwise) at the site j, 

3 

is the number operator for a spin-projection electrons, and t, U, and /zo are the first-neighbor transfer integral, the 
onsitc Coulomb repulsion, and the Bohr magneton, respectively. 

There are iVf up-spin electrons and N± down-spin electrons in the chain of N a sites and with lattice constant a 
associated with the model (1). We use periodic boundary conditions and consider iV a to be even and when N = N a 
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(half filling) both iVj and N± to be odd and employ units such that a = t = fiQ = h = l. The dimensionless onsite 
repulsion is u — U /At. When N a is odd the U = Fermi momenta are given by kp a = ± 



where 



k Fa = 77 ■ ( 6 ) 



N a 

When N a is even the Fermi momenta are given by hp = kp a and kp a = — kp a — jf- or by kp a = 



u„ _ _£ZL 

"■Fa Na 



and k Fa = —kp a . Often we can ignore the corrections of these expressions and consider kp a ~ ikp a — ±nn a 
and kp — [fe^-p + = irn/2, where 7v = N a /N a and n = N/N a - We consider the electronic density n = nj + nj. 

and the spin density m = n-f — rij,. 

For finite values of both the magnetic field and chemical potential the symmetry of the quantum problem is reduced 
to [7(1) (g) C7(l), with 5 C Z and S s z commuting with H |20| , |36| . The eigenvalues S c z and S s z determine the different 
symmetries of the Hamiltonian (1). When S c z ^ and S s z ^ the symmetry is U(l) ® C/(l), for = and 7^ 
(and = 0) it is 5?7(2) <g> C/(l), when S£ ^ and S* = it is [7(1) ® 5C/(2), and at S c z = and Sf = (and ^ = 0) 
the Hamiltonian symmetry is SO(A) . 

There are four U(l) <g> U(l) sectors of parameter space which correspond to S z < and S z < 0; S z < and S z > 0; 
S z > and 5f < 0; and S z > and 5f > 0. We follow Ref. Q and refer these sectors by (Z c , / s ) where 

Da 

■ ' (7) 
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The sectors (—1,-1); (—1,1); (1,-1); and (1,1) refer to electronic densities and spin densities < n < 1 and 
0<m<n;0<n<l and -n < m < 0; 1 < n < 2 and < m < (2 - n)\ and 1 < n < 2 and -(2 - n) < m < 0, 
respectively. 

There are two (l s ) sectors of SU (2)<g)[/(l) Hamiltonian symmetry [and two (l c ) sectors of U(1)®SU (2) Hamiltonian 
symmetry] which correspond to 5| < and 5f > for i s = — 1 and Z s = 1, respectively, (and to S% < and S 1 ^ > 
for l c — — 1 and Z c = 1, respectively). There is one 5*0(4) sector of parameter space [which is constituted only by the 
S z = (and [i = 0) and S z = canonical ensemble] . 

The related holon and spinon study of Ref. || did not introduce the excitation operator generators (some of these 
were studied in Ref. |20|] ) and only refers and applies to states with eigenvalues jf- ~ and ^f- ~ 0. In this paper we 
generate all Hamiltonian eigenstates from the pseudohole vacuum \V) introduced in Ref. p0|| , which is the half- filled 
and zero-magnetic-ficld GS. It is convenient for our complete operator description of the quantum problem to consider 
explicitly all the above sectors of parameter space. 

In the case of integrable models of simple Abelian U(l) symmetry the elementary excitations are generated by a 
single type of pseudoparticles (pseudoholes) [0. On the other hand, in the present case of the Hubbard chain we 
have shown Q that in terms of pseudoholes the description of the states I involves four branches of a, (3 pseudoholes, 
where a — c,s and (3 = ±i. In the particular case of LWS's I and (or) HWS's I we have that the a, (3 pseudoholes 
are such that = Therefore, in each (l c , l s ) sector of parameter space of Hamiltonian symmetry U(l) ® U(l) 
only the c, ^ and s, h£ pseudohole branches are involved in the description of the corresponding states I. However, the 
relation /3 = ^ is not valid for the general case considered in this paper, as we find below, yet the general description 
also involves the a, (3 pseudoholes introduced in Ref. pofl , which have creation operators a q a a. In addition, it involves 

an infinite number of a, 7 heavy-pseudoparticle branches, with creation operators 6j a „. Here the quantum numbers 
7 are the positive integers (7 = 1,2, ...,00). 

Elsewhere we will characterize the connection of the pseudohole and heavy-pseudparticle quantum numbers to 
electronic quantities in some limiting cases of physical interest. For instance, in the n — * 1 and S c — » (or m — » 
and 5 s — ► 0) limit, a — c (or a — s) becomes charge (or spin) and the c,/3 (or s, /3) pseudoholes are such that 
(3 = S z = ±| (or /3 = — 5*| = ^|) (20|. [We emphasize that in that limit and in the particular case we just considered 
5° = ±| (or 5| = ±1) denotes the 77 spin (or spin) projection of the individual c, f3 (or s, (3) pseudohole. Therefore, it 
does not refer to the general definition associated with Eqs. (3) and (4) which defines the rj spin (or spin) projection 
of an Hamiltonian eigenstate.] Moreover, for U — > 00, the creation of one c, 7 pseudoparticle (7 > 0) is associated 
with the creation of D = 7 doubly occupied sites. Note that in that limit the Hamiltonian commutes with the 
doubly-occupancy operator, i.e. [H, D] = 0. 

The above operators obey an anticommuting algebra such that 

{ a q,a,f3' a q',a',P'} = Sq,q' 5 a ,a' , (8) 
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and 



(9) 



and all remaining anticommutators vanish. 

As was found in Ref. j20|], the a, /3_rjscudoholes are the "holes" of the a pseudoparticlcs considered in the low-energy 

Till and studied in more detail in Refs. |32|,[33 34 3^]3(|. In contrast to their 



representation introduced in Refs 

"holes" , the usual a pseudoparticles (of number N a | 
p0|] . Moreover, there are no (3 bands, i.e. a, +5 and a 

the same a band [as the right-hand side (rhs) of Eq. (23) below confirms]. For convenience we shall call the usual a 
pseudoparticles p2 33 13~|,|35||3(| and the corresponding a bands the a, pseudoparticles and a,0 bands, respectively. 



56[) are not labeled by the (3 pseudohole quantum number 
i pseudoholes occupy different pseudomomentum values in 



The description of the whole Hilbert space involves the set of a, 7 pseudomomentum bands which are occupied by 
the a, (3 pseudoholes (or a, pseudoparticles) for 7 = and by the a, 7 heavy pseudoparticles for 7 > 0. In particular, 
the description of the non-LWS's and non-HWS's out of the BA requires the use of the a, (3 pseudoholes instead of 
the usual a, pseudoparticles. However, we find in futures sections that both the SO(A) Hamiltonian term (2) and 
the BA rapidities are (3 independent. Moreover, we find that the BA rapidities can be generalized to describe a whole 
SU(2) tower of r\ spin or spin Hamiltonian eigenstates. We emphasize that the BA solution refers directly only to 
the LWS's (or HWS's) of these towers p3fl . In the case of such /3-independent quantities it is often convenient to use 
the a, O-pseudoparticle representation rather than the a, /3-pseudohole representation. The a,0 pseudoparticles have 
creation operators b a and obey an anticommuting algebra such that 



(10) 



and all remaining anticommutators vanish (including the ones involving a,0 pseudoparticles and heavy pseudoparti- 
cles). 

Often, when we refer a, 7 pseudoparticles (bands) we consider 7 = 0,1,2,... whereas the notation a, 7 heavy 
pseudoparticles refers only to 7 = 1,2,.... The use of a, (3 pseudoholes and a,0 pseudoparticles is alternative and 
in general we do not use them at the same time. While in the characterization of all Hamiltonian eigenstates and 
in studies of /3-dependent operators the use of the a, /3-pseudohole and a, 7-heavy-pseudoparticle representation is 
required, in studies involving /3-independent operators (as for instance in charge transport |25|) it is often more 
convenient to use the a, 7-pseudoparticle representation (with 7 = 0, 1, 2, 3...). 

A crucial point is that the pseudohole and heavy-pseudoparticle numbers are constant for each Hamiltonian eigen- 
statc. The study of the non-LWS's outside the BA solution |23 requires the introduction of the generalized a, (3 
pseudohole numbers. These are given by 



N c,fi 



(3[N a - N] 



s,0 



/3[iV T - N{\ 



(11) 



where 



Nj} = N a ~ N c , , 



Aj 1 = A c , - 2 



7=0 



N, 



s,7 5 



(12) 



and the total number of a pseudoholes is N£ — . 

These important expressions are generalizations of the pseudohole expressions of Ref. |20j| , which refer only to the 
S c and S s LWS's and HWS's I (which we call simply states I). It is convenient to introduce the numbers N* such 
that 



N* = S a - l a S: = S a - \S?\ , (13) 

where for S" ^ the parameter l a = ±1 is given in Eq. (7). When dealing with energies and eigenvalues instead of 
Hamiltonians and operators the numbers of Eq. (13) together with the a, 7-pseudoparticle quantities and numbers 
constitute an alternative representation for the a, /3-pseudohole and a, 7-heavy-pseudoparticle quantities and numbers 
p6[ . The knowledge of all the N^ „ pseudohole and A Qi7 heavy-pseudoparticle (7 > 0) numbers provides the knowledge 
of all pseudoparticle numbers N an (7 = 0,1,2,3...) and numbers, and the inverse is also true. While the set of 
N an numbers are directly given by the BA solution Jl9[ , the /3-dependent pseudohole numbers ^ and the numbers 
are not considered by that solution. In Appendix A we relate the pseudoparticle numbers N an with the usual BA 
notation p9| . On the other hand, we find below that (3 is the quantum number needed to describe the Hamiltonian 
eigenstates out of the BA solution. 
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It follows from Eq. (11) that the electron numbers are exclusive functions of the pseudohole numbers and read 

Ni ; = ^y+J2 ^ - N U > N i = *y + E Pi**? + N "J ■ ( 14 ) 



As for the case of the Hilbert subspace spanned by the states I, the pseudomomentum-number operators 

N^p(q) = a\ a ^a q ^ p , N a „(q) = b\ arf b q , an , (15) 

(with 7 = 1,2, ...), and 

4o(?) = l-E^j(9), (16) 



play a central role in the generalized theory. The operator (16) has the following alternative representation in terms 
of a, pseudoparticle operators 

The number operators can be expressed in terms of the pseudomomcntum distributions as follows 

= E > Ki = E ^(<?) ■ as) 

We emphasize that while all Hamiltonian eigenstates are also eigenstates of the operator N£ ^ because these numbers 

are good quantum numbers, this does not hold true for the operator N% Jq). We find below that only the states I are 
eigenstates of the latter operator, yet all Hamiltonian eigenstates are also eigenstates of the pseudoparticle operators 

The 50(4) operator generators, S" and S± [see Eqs. (3) — (4)], and the ?y-spin and spin numbers S a and N* [see 
Eq. (13)], are given by 

s? = X>#£,/»(9); £2 = E a L±i c Wi> (1 9 ) 

q,0 1 

and 

1 V X 

S a = 2 [N* - £ ^N an ] ; N* a = -[^(1 - 2(3l a )N^ p - £ 2 7 iV Q , 7 ] , (20) 

7=1 (i 7=1 

respectively. Note that the expressions for the off-diagonal 5*0(4) generators go beyond BA, which provides the 
expressions only for diagonal operators. The application a suitable number of times of these operators onto the 
LWS's of the fy-spin and spin algebras generates the non-LWS's towers of Hamiltonian eigenstates outside the BA 
p3| . Therefore, their off-diagonal form in the quantum number /3, Eq. (19), reveals that in the present operator basis 
the construction of the states out of the BA requires the introduction of these pseudohole quantum numbers. This is 
confirmed by the form of the Hamiltonian-eigenstate pseudohole and heavy-pseudoparticle generators relative to the 
theory vacuum [see Eq. (23) below]. 

The generalization of the LWS's BA total momentum expression ||,[l9|], which we denote by Pba, to all Hamiltonian 
eigenstates, requires addition of an extra term associated with rj pairing |22|,^3|]. This leads to 

P = Pba + vr[5 c + S c z ] . (21) 

When this expression provides values of the momentum such that \P\ > ir, the total momentum is defined as the 
corresponding value at the first Brillouin zone. In operator form and for sub-canonical ensembles such that Sf 

and with l a given by Eq. (7) we have that N h _ la < N h la and the present representation leads to the following 

a, 2 a i — 

simple expressions for the momentum operator P and its eigenvalue P 

00 

P = J2 E Q^N^iq) + n[N^ + E ; 

g,CK 7— 7=1 
OO 

P = E E iC^N^q) + 7T [(1 + 7)^ c ,7 + K } , (22) 

q,a 7 — O 7 — 1 



G 



respectively, where C c . 7 = —1 for 7 > and C an — 1 otherwise. We emphasize that the momentum term, 
7r S 7 =i [(1 + 7)^c,7 + N£], is always a multiple of ±ir. 

The pseudoparticle perturbation theory introduced in Refs. 0,|],|]j3|]3^||||34| and developed in a suitable operator 
basis in Refs. ]s|j20|,|36|l refers to the Hilbert subspace spanned by the Hamiltonian eigenstates I. Rather than holons 
and spinons p|, |20| , |27| , at finite values of the magnetic field and chemical potential and at constant electronic numbers 
the low-energy excitations I are described by pseudoparticle-pseudohole processes relative to the canonical-ensemble 
GS. The latter state, as well as all excited states I with the same electron numbers, are simple Slater determinants of 
pseudoparticle levels P, |2fj| , p8| , p6[ . 

Although we find in Sec. Ill that the pseudoholes and heavy pseudoparticles obey generalized ^| rather than Fermi 
or Bose statistics, one can use their formal anticommutation relations (together with suitable topological-momentum- 
shift operators introduced in Refs. [^0 30 and generalized for the heavy pseudoparticles in Sec. Ill) to develop an 
extremely simple second-quantized representation of all 4 Na orthonormal Hamiltonian eigenstates |23|| . As for the 
above states I, in our generalized operator basis all these states are simple Slater determinants of pseudohole and 
heavy-pseudoparticle levels. They are of the form 



1 00 
\t,{N^},{N a , 7 }) = ^l[A a l[ 



[SI] 



9,9' 



■i q',a,~f 



\V), 



where the generator is given in Eq. (19), {N% p} and {N an } are abbreviations for the sets 

l., 2 2 a, 2 o, 2 

and 



(23) 



(24) 



{N Cjl ,...,N ( 



N 1 

c,oo ? 1 v s,l ; 



N \ 

) 1 v S,OC / 1 



respectively, 



C 



n( 



1 Y a • 

N h , 



(3 



7 = 1 



(25) 



(26) 



and 6(x) = 1 for x > and Q(x) = for x < 0. Note that A Q equals either 1 or 0. This is because in the present 
representation the Hamiltonian eigenstates are characterized by pseudohole and heavy-pseudoparticle numbers such 
that N£ p > Yl^Li 7-^a,7- We note that application of the non-diagonal opertors of Eq. (19) onto i] spin or spin 
LWS's or HWS's generates often states with numbers such that N% p < E^i7^a,7- The presence of A a (see Eq. 
(26)) in the rhs of Eq. (23) (which is to be changed according to the changes of the numbers Nj£ p) then assures that 
application of the above generators in these LWS's or HWS's gives 0. 

The existence of the half-filling Mott-Hubbard gap [^8|,^2| , which we denote here by Amh , implies that S c — — 
for chemical potentials fj, such that - Amh. < M < Amjl_ Therefore, the S c = S c z = and S s = S s z = vacuum \V) of 



the rhs of Eq. (23) can be associated with different chemical-potential values within the range 



< n < 



2 1 > ■ 

The vacuum chemical-potential value is choosen according to the suitable reference level from which the energy of 
the pseudohole and pseudoparticle seas of the states (23) is to be measured. Therefore, we assume in the definition 

. Ajv j w r — S „ 



of the vacuum of Eq. (23) that \V) is such that /1 = for S c z — states but /u — ► ±- 



for 



The possible pseudohole and heavy-pseudoparticle numbers are constrained by Eqs. (19) and (20) and together 
with the number of discrete pseudomomentum values in each band leads to A Na possible orthonormal Hamiltonian 
eigenstates of form (23), in agreement with Ref. [ p3[ . The number of discrete pseudomomentum values in each a, 7 
band (a, 7-band Fock-space dimension) is given by 



d Fan =N>: + N an - Y, [7 + 7' - |7 - 7'|]A^, 7 

7 ' = 1 

In the above summations and products, the pseudomomentum takes on values 

2t 7-0.7 



(27) 



(28) 



7 



where in contrast to the usual momentum, J"' 7 are consecutive integers or half-odd integers for iV Qi7 odd and even, 
respectively. Here 

N h °° N 00 

^c,o = E^-E^7 = ^-^,o- E n °,t ( 29 ) 



7=1 0,7 = 1 



and 



iV Q , 7 = d Fa , 7 , (30) 



or 



for values of a, 7 other than c, 0. It follows that for each a, 7 band < q < tta/yj with the limits of the pseudo- 
Brillouin zones given by 

<&=±^~^fh (31) 

for iV C! o even and 

8$=*, q$ = -*[i-^- a ], (32) 

?S } =T[l-j^] 5 (33) 
for iV Ci o odd. On the other hand, the limits of the pseudo-Brillouin zones are simply given by 

q&=±q a „ = ±-£r[dF a , y -l], (34) 

for all the remaining a, 7 bands (other than the c, band). 

The usual first-quantized BA wave function expression for the Hamiltonian eigenstates has an involved form which 
includes many permutations jl8|,[l9]]. That spatial wave function depends on the quantum numbers J"' 7 of Eq. (28) 
through infinite sets of complex numbers, which many authors call rapidities. The expression of the spatial wave 
functions in terms of the quantum numbers J"' 7 requires the solution of an infinite system of algebraic equations 
which we present in Sec. IV. These equations define the real part of the BA rapidities as functions of the quantum 
numbers If 1 @|f]. Although the expression of the spatial wave function for the Hamiltonian eigenstates in terms 
of the quantum numbers J"' 7 requires the solution of the above systems of equations, which constitutes a problem of 
considerable complexity, the description of these eigenstates in the basis of the above BA quantum numbers I"' 7 that 
diagonalize the quantum liquid, Eq. (23), is much simpler. Historically, the Hamiltonian eigenstates were introduced 
in terms of the spatial BA wave functions Jl||,[l5|] . The diagonalization of the problem leads then to infinite systems 
of algebraic equations. These equations introduce the integer or half-odd integer quantum numbers I"' 1 of Eq. (28) 
which describe the Hamiltonian eigenstates. 

One of the principal advantages of the representation (23) is that it permits the description of the Hamiltonian 
eigenstates in terms of the quantum numbers J"' 7 and does not require the spatial wave-function representation. In 
the basis associated with these quantum numbers the description of these states does not involve the rapidity numbers. 
Moreover, our algebraic representation (23) refers also to the Hamiltonian eigenstates out of the BA solution. 

A S a LWS is such that 



whereas for a S a HWS we have that 



<±=E^7, (35) 
7=1 



7=1 



with N h ! > N h 1 for a LWS and N h 1 > N h t for a HWS. 
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For the particular case of states I, we have that iV a . 7 = for all 7 = 1,2, ... branches in the former equations. 
Therefore, the states I are characterized by zero occupancies of the heavy-pseudoparticle bands a, 7 (with 7 > 0). It 
then follows from the general Eqs. (35) and (36) that a S a LWS I is such that 

^i=0, (37) 

whereas for a S a HWS I we have that 

tf£_i=0, (38) 

' 2 

with A^ l _i > for a LWS I and t > for a HWS I. This together with Eq. (23) leads to the simple Slater- 
determinant form found for the states I in Ref. p0|| . These results also confirm that in the case of states I the a, f3 
pseudoholes are such that /3 = but that this equality does not hold true in the general case. While in the case of 
the states I the Slater determinant (23) either involves a, +4 or a, —4 pseudoholes only, in the case of the non-LWS's 
and non-HWS's or in the case of states II (i.e. LWS's or HWS's containing a, 7 heavy pseudoparticles) the state (23) 
has both cx,+^ and a,— 4 pseudohole occupancies. 

In the case of the non-LWS's and non-HWS's out of BA the relative number of a, \ and 0^,-5 pseudoholes is 
determined by the values of S a and S" through Eqs. (19), (20), and (39). The pseudohole numbers which characterize 
these states are such that 

00 

- £ 2 7 iV Q , 7 > \N* ti - N*_ 1 1 . (39) 

7=1 

They can be generated from LWS's (or HWS's) by acting onto the latter the off-diagonal generator S" (or S") the 
required number of times. According to Eq. (19), in the present basis these off-diagonal generators just produce /3-flip 
processes. In the case of the LWS's II the (3 flips generated by these operators create a number 7 of a, ^ pseudoholes 
for each a, 7 heavy pseudoparticle, which implies these states are inside BA. Following Eqs. (19) — (20) this is required 
for the LWS condition S a = -S?. Note that both for the non-LWS's and non-HWS's out of the BA and for the 
LWS's II and HWS's II there is no relation between the quantity l a defined in Eq. (7) and the pseudohole quantum 
number (3. 

In the case of non-LWS's and non-HWS's the wave function (23) uses as starting state a LWS. We then call it for 
LWS general representation. The general state (23) can be rewritten, alternatively, in terms of a starting HWS (HWS 
general representation). This just requires replacing in the rhs of Eq. (23) ±4 by and the generator S" by 5". In 
the corresponding HWS general representation the off-diagonal generators create new states from a reference HWS. 
One can also use a representation where the starting state is a S c (or S s ) HWS and a S s (or S c ) LHS. However, the 
four possible general representations describe the same states. 

The form of the pseudohole and heavy-pseudoparticle generators of the rhs of Eq. (23) confirms that the correspond- 
ing Hamiltonian eigenstates are not in general eigenstates of the pseudohole operator Jq). This is because due to 
the transformation laws of the off-diagonal fy-spin and spin generators present in the rhs of Eq. (23) a Hamiltonian 
eigenstate is, in general, a mixture of different (3 — — 4 and j3 — +5 pseudomomentum distributions. The exception 
are the states I which have only either a, \ 01 q,-| pseudohole occupancy. 



III. GENERALIZED GROUND STATES AND STATISTICS 



The conservation of the pseudohole and heavy-pseudoparticle numbers permits dividing the Hilbert space into 
subspaces spanned by the set of Hamiltonian eigenstates (23) with the same g} and {N a>1 } numbers, which 
correspond to the same sub-canonical ensemble. Obviously, a canonical ensemble (with constant electron numbers 
and thus with constant values of N 1 ^ L — _ L for both a = c,s) is usually realized by several {N^A, {A^ Q . 7 } 

sub-canonical ensembles. Let us introduce the generalized ground state (GGS) as the Hamiltonian eigenstate(s) (23) 
of lowest energy in each Hilbert subpace. We find it is (they are) of the form 



\GGS;{N^},{N an }) 



1 OO 

v a 7—1 



* a + f -* n n n n n 

9=9i, ?=9fc,,o 9'=9c,7 9'=9f^! 7 l"~lFB, n 



9f q ,0 



9Fc,7 



V,-r« , ,=,7T 1 «,7 



\v) 



(40) 



9 



where 



(±) (±) ±r 2vr 

and when (i) a = c, 7 > 0, and -/V Cj7 is even or (ii) a = s or a = c and 7 = and N a ,-y is odd (even) and 7"' 7 are 
integers (half-odd integers) the pseudo-Fermi points are symmetric and given by 

q l F^=±[qF a ,y-C Ch ^], (42) 



where 



?Fc, 7 = ^7 ; 7 > U , 



QFc,0 



N a 



<?F S , 7 = TT ■ (43) 

On the other hand, when (i) a = c, 7 > 0, and AT c-7 is odd or (ii) a = s or a = c and 7 = and N an is odd (even) 
and I"' 7 are half-odd integers (integers) we have that 

iF^.j = lFa,j , 9Fi 7 = _ [9F Q ,7 - C Qi7 — ] , (44) 



ffF^ = 9Fa,7 - C a , 7 — , g^J 7 = -gFa, 7 • (45) 



The GS associated with a given canonical ensemble is always a state I of the form presented in Ref. J20|, which is a 
particular case of the general GGS expression (40). It is useful to denote the GS's by \GS; S%, S%), where S% and 5| 
are the eigenvalues of the corresponding canonical ensemble. The GS expression associated with the (l c , l s ) sector of 
Hamiltonian symmetry U(l) ® U(l) reads 



=<-) <+) =(-) (+) 



\Gs;st,si)= n n < c , ¥ n n i, s ,^- ^ 

1=li~o 1=<1 { fIo 1=1^0 1=^,0 

In all sectors of Hamiltonian symmetry there are states I. In the particular case of the SO (4) zero-chemical potential 
and zero-magnetic field canonical ensemble there is only one state I, which is the pseudohole and heavy-pseudoparticle 
vacuum, \V), which is nothing but the SO(4) GS [P0[ . The same applies to the sectors of Hamiltonian symmetry 
SU(2) ® U(l) and U(l) (8> SU(2), the GS being always a state I. (In addition, in these sectors there is a large number 
of excited states I.) 

In the case of the two (l a ) SU(2) ® U(l) sectors the GS is both a LWS and HWS of the 77-spin algebra. Therefore, 
it is empty of c pseudoholes and reads 

\gs;o,s° z )= f[ n i^jvy ( 4? ) 

In the case of the (l c ) U(l) ® SU(2) sector the GS is both a LWS and a HWS of the spin algebra and is empty of 
s pseudoholes. It reads 

|GS;S 2 c ,0) = [] II a Ic^l y )- ( 48 ) 
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Finally, the S c = = (and /i = 0) and S" = S* = 50(4) ground state is, at the same time, a LWS and 
HWS of both the 77-spin and spin algebras, ie following the notation of Refs. @|| it is a [LWS,LWS], a [LWS,HWS], 
a [HWS, LWS], and a [HWS, HWS]. Therefore, it is empty of both c and s pseudoholes and is the vacuum of the 
pseudohole and heavy-pseudoparticle theory. All the remaining states I can be described by particular cases of the 
general expression (23) which are Slater determinants of a, (3 pseudoholes and a, 7 heavy pseudoparticles pQ] , 

The number of Hamiltonian cigcnstates associated with a sub-canonical ensemble and thus characterized by the 



same {JV'J and {N an } numbers is given by 



n {iV a , 7 }) = J] 



d-FaA 



a,+4 



N 



n 

7=1 



d-Fo 



n 



dFafi 

N h 



1 00 r 

n 

7=1 



which in terms of pseudoparticle numbers only can be rewritten simply as 



n n 

a— c,s 7— 



(49) 



(50) 



The square brackets in the above equations refer to the usual combinatoric coefficients. 

To classify the statistics of our pseudoholes and heavy pseudoparticles according to the generalized Pauli principle 
of Ref . fl29j , we introduce the pseudohole and heavy-pseudoparticle dimensions 



1 



N h 



1 



(51) 



respectively. An Hamiltonian eigenstate transition producing changes AJV* and AiV Q . 7 in the pseudohole and heavy 
pseudoparticle numbers, respectively, leads then to the following changes in the corresponding dimensions 



a' a' 7' — ! 



(52) 



and 



a.' a' 7' — 1 

where the statistical-interaction matrix has infinite dimension and reads 

9c,oc — ^c,q: 3 Pc;a,7 ) 



(53) 



(54) 



1 

2 ' 



^s;a,7 — $s,a j 



(55) 



This fully defines the statistics, the pseudoholes and heavy pseudoparticles being neither fermions nor bosons |29| ]. 

The simple form of the GGS expression (40), GS expressions (46) — (48), and general-Hamiltonian eigenstates (23) 
has a deep physical meaning. It reveals that in the present basis these eigenstates of the many-electron quantum 
problem are "non-interacting" states of simple Slater-determinant form. However, that the numbers I"' 7 of the rhs 
of Eq. (28) can be integers or half-odd integers for different Hamiltonian eigenstates, makes the problem much more 
involved than a simple non-interacting case. This change in the integer or half-odd integer character of some of the 
numbers I"' 7 of two states, shifts all the occupied pseudomomenta and leads to the orthogonal catastrophes ||. These 
shifts are associated with topological excitations which are generated by the topolological-momentum shift operators 
introduced in Refs. fl20| , |3"cfl for the particular case of the states I. 

Let us introduce the a,0 and a, 7 topological-momentum-shift operators (for the states I see Refs. po| , [30f ) which 
generate topological excitations associated with a shift of all the pseudomomenta of the corresponding band by ±t5-. 
Creation or annihilation of odd numbers of pseudoholes or heavy-pseudoparticles always requires the occurrence of 
such topological-momentum shifts. Therefore, in addition to the pseudohole and heavy-pseudoparticle operators our 
second-quantized algebra also includes topological-momentum shift operators. Their unitary generators are such that 
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^0' 



±1 



(57) 



and 



TT±l h t =h t TT±1 



(58) 



and read 



<7,/3 



(59) 



and 



(60) 



As we discuss in future sections, a remarkable property is that the transition between a GS (46) and any eigenstate 
(23) can be separated into two types of excitations: (a) a topological GS - GGS transition which involves the creation 
or annihilation of pseudoholes and (or) heavy pseudoparticles as well as the occurrence of topological momentum 
shifts and (b) a Landau-liquid excitation ||,^0| associated with pseudoparticle - pseudohole processes relative to the 
GGS. 

The topological transitions (a) are basically superpositions of three kinds of elementary transitions: (i) GS - GS 
transitions involving changes in the a electron numbers by one, (ii) single (3 flip processes of a, (3 pseudoholes which 
lead to non-LWS's and non-HWS's outside the BA, and (hi) creation of single a, 7 pseudoparticles at constant values 
of S c , S^, S s , and 5|. While the transitions (i) are gapless, the elementary excitations (ii) and (hi) require a hnitc 
amount of energy. 

The generators of the excitations (i) were studied in Ref . p0| and the ones of (ii) involve only topological momentum 
shifts and (3 flips (which describe either creation of electron pairs or spin flip processes). Consider as an example 

of a transition (hi) the creation of one a, 7 pseudoparticle at the pseudo- Fermi points, q = q^^- (We examine 
here explicitly topological momentum shifts of occupied bands only, i.e. these which generate momentum.) Let us 
introduce the notation —a such that —c = s and — s = c. When 7 > is even (or odd) the generator reads U^\G atl 
(or G a ~U^a ), where the topological-momentum-shift operator was defined above and 

(-) + 2 JL , ^ ( + ) 

G a „=[str n n a L^( ± ), Q , 7 - ( 6i ) 

a=a ( ~ ) o <+) -^-(-y-1) 

Therefore, such transitions involve one a, or —a, topological momentum shift, the creation of a number 7 of a, h 
pseudoholes and 7 of a, — ^ pseudoholes, and the creation of one a, 7 pseudoparticle at one of the pseudo- Fermi points. 
(This can be considered as the antiparticle of the "cloud" of 27 associated pseudoholes.) 



IV. PSEUDOPARTICLE HAMILTONIAN AND BA RAPIDITIES 



For simplicity let us denote the general Hamiltonian eigenstates {N^ g}, {N an }) by 1^) [and the GGS's (40) by 
\GGS).] These states arc also eigenstates of the a, 7 pseudomomentum-distribution operators (15) such that 

IV) = N a „(q)\il>) , (62) 

where N a ,j(q) represents the eigenvalue of the operator (15), which is given by 1 and for occupied and empty values 
of respectively. 

It follows from the form of the state-generator of the rhs of Eq. (23) that all 2S a + 1 Hamiltonian eigenstates 
constructed from the same S a = —S" LWS by applying onto it 1, 2, 2S a number of times the off-diagonal operator 
S" (19) have the same N an (q) pseudomomentum distribution. On the other hand, the pseudohole numbers _i 

and N h , 1 are different for each of the 2S a + 1 states of such S a tower. These numbers change from TV' 1 1 = N% — 
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E~= i 7^a, 7 and N h a +h = £~ =1 7 iV Q , 7 for the initial LWS to N^_, = £" =1 7^,7 and 2V£ +J = - E"=i 7^,7 

for the corresponding HWS, where iV^ is given in Eq. (12). 

As in the case of the spatial wave functions, the expression of the energy in terms of the quantum numbers I"' 7 of Eq. 
(28) involves the BA rapidities. These are complex functions of the numbers ^"' 7 - Following our pseudomomentum 
definition, Eq. (28), in our description they arc functions of the pseudomomentum q. The real part of the BA rapidities, 
which we denote by i? Q;7 (g), are in the present basis eigenvalues of corresponding rapidity operators, i? Qi7 (q). If we 
combine our operator representation with the properties of the BA solution, it is straightforward to confirm from 
the relation between the rapidity real part, R an (q), and the pseudomomentum distributions of Eq. (62) that these 
real rapidity- function parts are eigenvalues of the rapidity operators R a (q). These rapidity operators are exclusive 
functions of the pseudoparticle operators N an (q). It follows that all Hamiltonian eigenstates are also eigenstates of 
the rapidity operators R a ^(q) such that 

R a , 7 (q)\iP) = R a , 7 (q)\-il>) , (63) 

where for the particular case of R C fi(q) we also consider the associate rapidity K(q) which is defined by the following 
equation 

„ , , sin K(q) 

R c .o(q) = — • (64) 

u 

This notation refers to the spectral parameter ESjHgi which appears often in operator form. The rapidity K(q) is 
the eigenvalue of the rapidity operator K(q) such that 

K(qM = K(q)\i/>). (65) 

In Appendix A we present the relation of our rapidity eigenvalues to the real part of the B A rapidity numbers [fl9[ . 

The real part of the BA rapidities, which in our basis are the eigenvalues R an (q) [and K(q)], are fully determined 
by the pseudoparticle distributions N a ^(q) through a system of coupled integral equations. The standard treatment 
of the BA solutions is very lengthy, e.g., for each state we have to rewrite a new set of equations. This is because the 
eigenvalues i? aj7 (q) of the rapidity operators R a:7 (q) are different for each state \ip) (23). 

However, within the present operator description, we can introduce a single set of infinite general equations which 
apply to all LWS's. In addition, our operator representation allows the generalization of the BA rapidities to the 
non-LWS's out of BA. It follows that the above general equations can be generalized to all Hamiltonian eigenstates. 
Each of these equations defines one of the rapidity operators R an (q) in terms of a q, a, and 7 summation containing 
functionals of the rapidity operators and pseudomomentum-distribution operators. This system of coupled equations 
has a unique solution which defines the expressions for the rapidity operators in terms of the pseudomomentum- 
distribution operators. In terms of the corresponding eigenvalues the system of coupled integral equation reads 

(<Z) ^ _ A^^ s ' 7(?) \ 1 + 7 

° 7=0 q' 

' W E E *o,M) tan- , (66) 

7=1 q> 

2Re sin" 1 (u(R c , 7 (q) - i 7 )) = q + J- £ N c , (q>) tan" 1 ~ Rc ' o{q>) ' 

a <?' 

+ w EE^'feOe^'t^W - Rc,y(q')) , (67) 

° 7' = 1 q' 

and 

= q - J- £ N M) tan" 1 ( p^Mtl^M 

a q> ' 



r EE^^Hi+ykta)-^^)) . (68) 



7'=0 q' 
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where 



7-1 

©7,7' (*) = ©7',7 (*) = ^7,7'{2 tan" 1 (|-) + ^ 4 tan" 1 (J) } 



Z=l 



,2 7 

+ (1 - * 7 ,,){2 tan- 1 ( ] _£- ;[ ) +2 tan" 1 (-JL_) 

2+7 — I 

+ | "--'(j^pi+s)}- (») 

The limits of the pseudo-Brillouin zones, q^y, associated with the pseudomomentum summations are given in Eqs. 
(31) — (34). Although the integral Eqs. (66) — (68) are coupled note that each of these equations refers to one a, 7 
channel and defines the corresponding a, 7 rapidity in terms of other rapidities. 

The key point is that since the eigenvalue iV Qi7 (g) (here 7 — 0, 1, 2, 3, ...) is common to the whole tower of 2S a + 1 
Hamiltonian eigenstates constructed from the same S a LWS, also the associate rapidity eigenvalue i? Qj7 (g) of Eq. 
(63) is common to these 2S a + 1 Hamiltonian eigenstates. This follows from the fact that rapidity solutions R an (q) 
of Eqs. (66) — (68) are functionals of the N an (q) distributions and (3 independent. We thus conclude that the LWS 
rapidities extracted from the BA provide full information on the non-LWS's rapidity operators. 

The rapidity eigenvalues of Eqs. (63) and (65) are independent of the pseudohole numbers ^ and only involve the 
pseudoparticle distributions N an {q). The same holds true for the associate rapidity, pseudohole, and pseudoparticle 
operators. Since the 5*0(4) Hamiltonian (2) and other physical quantities turn out to be exclusive functionals of 
the rapidity operators, in studies involving such quantities it is often more convenient to use the a, 7 pseudoparticle 
representation (with 7 = 0,1,2,3,...) than the a, (3 pseudohole and a, 7 heavy-pseudoparticle representation (with 
7 = 1,2,3, ...). We remind that the set of all ^ and N U:7 (with 7 = 1,2,3, ...) numbers of a given Hamiltonian 
eigenstate (23) have the same information as the corresponding set of all N* and N a i7 (with 7 = 0, 1, 2, 3, ...) numbers 
which characterize the same state. 

By direct insertion in each of Eqs. (66) — (68) of the corresponding gi^ pseudomomenta we find that the functions 
R ai7 (q) have for all eigenstates the following boundary values at the limits of the pseudo-Brillouin zones 

K{q$>) = ±v, i? c , (^ ± ) )-0, (70) 

for a = c and 7 = and 

fla,7(O = ±0 °> ( 71 ) 

for the remaining choices of the quantum numbers a, 7 other than c, 0. 

It is useful to consider the rapidity functions K^°\q) and Ra)y(q) which are the GGS eigenvalues (or GS eigenvalues 
- a GS is a particular case of a GGS) such that 

K(q)\GGS) = K ia \q)\GGS) , (72) 

and 

R a . 7 (q)\GGS) = RM(q)\GGS) . (73) 

These are defined by Eqs. (66) — (68) with the pseudomomentum distribution given by its GGS value which we denote 
by Na}/(q) and is such that 

N an (q)\GGS) = NM{q)\GGS) . (74) 

Following the GGS expression (40), Na}/(q) has the simple free Fermi- like form 

JV<™ (q) - 6 (C a , 7 [<4 + j, 7 -q]), < q < 

= 9 (C ar/ [q - g^g) , g<-) < q < . (75) 

For the particular case of a GS these distributions read 



14 



= ^(q-9f1o)> lL~0<1<0: 



N<®(q)=0, 7 >0. (76) 

By inserting in Eqs. (66) — (68) the GGS distribution (75) we define the GGS rapidities K^(q) and Ra)y(q) in 
terms of their inverse functions with the result 

K<°>(g) 

dk2irp c , (k) , (77) 



dr2Trp a ^{r), (78) 

where the fuctions 2np C fi{k) and 27rp Q . 7 (r) are the solutions of the integral equations (B4)-(B6) of Appendix B. In 
that Appendix we study the rapidity solutions of Eqs. (66) — (68) both for GGS's and for Hamiltonian eigenstates 
differing from a GGS by a small density of excited pseudoparticles. In equations (B4)-(B6) the parameters and 
Ta^y are defined combining the following equations 

Q {±) = ^ (0) (^o) , = R^Mf%) . (79) 

with equations (77) and (78). We also introduce the associate parameters 

c,o) i ^a,7 

= RW(q Fa!j ). (80) 

The system of coupled algebraic equations (66) — (68) is obtained from the BA solution (see Ref. |l!| and Appendix 
A which connects our rapidity representation to the notation used in that reference), yet also describes states out of 
that solution. Given the configuration of pseudoparticle quantum numbers that describes each Hamiltonian eigenstate 
of form (23), these equations fully define the corresponding rapidities. In the thermodynamic limit, we can take the 
pseudomomentum continuum limit qj — > q and the real part of the rapidities become functions of q which we have 
called here R an (q). In that limit the BA system of algebraic equations are replaced by the system of infinite coupled 
integral equations presented in Appendix B. 

We have mentioned that the BA solution is most naturally expressed in the pseudoparticle basis. One reflection of 
this is the simple expression for the Hamiltonian (1) in that basis. The BA energy expression given in Refs. |], |l9) for 
the LWS's of the 77-spin and spin algebras, which we denote by Eba, can be generalized to all Hamiltonian eigenstates. 
For that we replace the —5" dependence of Eba by a S a 

E = E BA (-S«^S a ). (81) 

Combining Eq. (19) with the energy expression (81) leads then to the following expression for the Hamiltonian (1) 

H = H so{4) + 2 M ]T 0^o,p + 2 Mo# E Pfis,p , (82) 

where 



H S o(i)=-2t^2N c , (q)cos[K(q)]+2t ]T N^(q) ]T yj 1 - [u(R c ^(q) ~ ji 7 )] 2 

q 9,7=1 i=±l 

+ ^-% £ -£^]- (83) 

7=1 

This is the rigurous pseudoparticle expression of the Hamiltonian (2) at all energy scales. It gives the exact expression 
of that Hamiltonian in the full Hilbert space. Despite its simple appearance, the Hamiltonian (83) describes a many- 
pseudoparticle problem. The reason is that the expression of the rapidity operator in terms of the operators N atl (q) 
contains many-pseudoparticle interacting terms. 
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To specify the Hamiltonian (83) completely, we must indicate the equations that define the rapidity operator 
eigenvalues in terms of the pseudoparticle momentum distribution operators. In the case of the present model these 
are Eqs. (66) - (68). 

Unsurprisingly, it is difficult to solve the BA operator equations (66) — (68) directly and to obtain the explicit 
expression for the rapidity operators in terms of the pseudomomentum distribution operators (15). In contrast, it is 
easier to calculate their normal-ordered expression in terms of the normal-ordered operators. The rapidity operators 
Ra,-y(q) contain all information about the many-pscudoparticle interactions of the quantum- liquid Hamiltonian. There 
are two fundamental properties which imply the central role that the rapidity operators of Eqs. (72) and (73) have 
in the present quantum problem: 

(a) As we find in Sec. V, each of the normal-ordered rapidity operator : R a (q) : (relative to the suitable GGS or 
GS) can be written exclusively in terms of the pseudomomentum distribution operators (15); 

(b) The normal-ordered version of the S'0(4) Hamiltonian (83) can be written, exclusively, in terms of the pseu- 
domomentum distribution operators (15), but all the corresponding many-pseudoparticle interaction terms can be 
written in terms of the rapidity operators : R a ^(q) :. It follows that the rapidity operators commute with the 
Hamiltonian. 

In the ensuing section, we introduce the pseudoparticle perturbation theory which leads to the normal-ordered 
expressions for the Hamiltonian and rapidity operators. Despite the non- interacting form of the Hamiltonian eigen- 
states (23), the normal-ordered Hamiltonian includes pseudoparticle interaction terms and is, therefore, a many- 
pseudoparticle operator. However, we find in Sec. V that these pseudoparticle interactions have a pure forward- 
scattering, zero-momentum transfer, character. 



V. PSEUDOPARTICLE HAMILTONIAN AND PERTURBATION THEORY 



We will be mostly interested in GS - GGS transitions followed by pseudoparticle - pseudohole excitations involving 
a small density of a, 7 pseudoparticles relative to both the GS and GGS distributions. Let us then introduce the 
normal-ordered pseudomomentum distribution operator 

■N an (q) := JV Qi7 (g) - JV<$ («) . (84) 
We emphasize that when we choose Najy{q) to be the GS pseudomomentum distribution (76) we have that 

: N a . ( {q) :=N an (q), 7 >0. (85) 

The normal-ordered distribution (84) obeys the eigenvalue equation 

: N a „(q) : \i/>) = 5N a „(q)\1,) . (86) 
We also consider the normal-ordered rapidity operator 

: R a , 7 (q) := R a , 7 (q) - R^(q) , (87) 

where R^ X J (q) is the GGS eigenvalue of Eq. (73). 

In the pseudoparticle basis the normal-ordered rapidity operators : R an (q) : contain an infinite number of terms, as 
we shall demonstrate below. The first of these terms is linear in the pseudomomentum distribution operator : N arf (q) :, 

(84), whereas the remaining terms consist of products of two, three, , until infinity, of these operators. The number 

of : N an (q) : operators which appears in these products equals the order of the scattering in the corresponding 
rapidity term. 

A remarkable property is that in the pseudoparticle basis the seemingly "non-perturbative" quantum liquids become 
perturbative: while the two-electron forward scattering amplitudes and vertices diverge, the two-pseudoparticle / 
functions [given by Eq. (117) below] are finite. 

By "perturbative" we also mean here the following: at each point of parameter space (canonical ensemble) the 
excited eigenstates arc of form (23) and correspond to quantum-number configurations involving a density of excited 
pseudoholes and heavy pseudoparticles relative to the GS configuration (46) and (76). We introduce the following 
density 

<* = E<^+ E ( 88 ) 

Q./3 a, 7=1 
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which is kept small. Here 

n a,/3 = 

N a 



<£ = ±- £[1 - N h a f(q)]SN a Aq) , (89) 



and 



,<*>7 
l ex 



define the densities of excited a, (3 pseudoholcs and a, 7 heavy pseudoparticles, respectively, associated with the 
Hamiltonian eigenstate \ip). 

In Eq. (89) A^'s (q) is the GS pseudohole distribution which reads 

N^l(q) = l-N^(q), 



N 



h,(0) 



= 0, (91) 



if the GS is a S a LWS and 



<i°](<?)=i-<<k?), (92) 

if the GS is a 5 Q HWS, where N^ (q) is the a, pseudoparticle GS distribution (76). Our perturbation theory refers 
to general GS transitions to Hamiltonian eigenstates involving (a) a GS - GGS transition which costs the energy 
u>o [given by Eq. (124) below] and (b) a Landau- liquid pseudoparticle-pseudohole excitation around the GGS. The 
pseudoparticle perturbation theory corresponds to such GS transitions with excitation energy lu such that the energy 
(lu — luq) is small. Since the initial state is a GS we have assumed the initial GS distribution (76) for the a, 7 heavy 
pseudoparticles. Following Eq. (76) we thus have that N a , 7 (q) = SN atl (q) in the rhs of Eq. (90). 

The expectation values of the S'0(4) Hamiltonian (83) in the final states \ip) are functions of the density of excited 
pseudoparticles only. This density is given by 



(93) 



a,7— 

where is given by 



n -~° = W a ^ [1 ~ JV S(9)]*JVa,o(?) , (94) 



<>.x 



for 7 = and by (90) otherwise. When both the pseudoparticle density (93) and n z = — - [see Eq. (13)] are 
small implies that the density (88) is also small. On the other hand, the pseudoparticle density (93) is small provided 
the density (88) is also small. This then implies that all the densities n"^ are small and that we can expand the 
expectation values in these densities. The perturbative character of the quantum liquid rests on the fact that the 
evaluation of the expectation values up to the j th order in the densities (93) requires considering only the corresponding 
operator terms of scattering orders less than or equal to j. This follows from the linearity of the density of excited 
a, 7 pseudoparticles, which are the elementary "particles" of the quantum liquid, in 

6N a „(q) = {tl>\iN a „(q):\1>), (95) 

and from the form of (90) and (94). The perturbative character of the quantum liquid implies, for example, that, to 
second order in the density of excited pseudoparticles, the energy involves only one- and two-pseudoparticle Hamil- 
tonian terms, as in the case of the quasiparticle terms of a Fermi-liquid energy functional In this case the 
corresponding energy (cj — uiq) is small. 

It follows from the eigenfunction Eq. (86) that the problem of using the rapidity Eqs. (66) — (68) to derive the 
expression of the operators : i? Qj7 ((?) : in terms of the operators : N an (q) : is equivalent to the problem of evaluating 
the corresponding expansion of the rapidity eigenvalues 5R an {q) in terms of the pseudoparticle deviations SN ail (q) 
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(95). This last problem, which leads to the Landau-liquid expansions, was studied in Refs. |3^|5j,|36| f° r the case of 
the Hilbert subspace spanned by the states I. Furthermore, we emphasize that it is the perturbative character of the 
pseudoparticle basis which justifies the validity of these Landau-liquid deviation expansions |20|]3^ , |3^ , |34] , |36t . 

Since a GS is a particular case of a GGS, in the following we consider the general case of normal ordering relative 
to a GGS. Only by using as reference state a GGS (with heavy-pseudoparticle occupancy) and obtaining the GS 
expressions as a particular limit of the GGS general expressions we can be shure to arrive to the correct GS results. 
Later on we will specify when the ordering refers to a GS. 

We can derive the expressions for the rapidity operators : R a ^(q) :. This corresponds to expanding the expressions 
of the operators : R a ^{q) : in terms of increasing pseudoparticle scattering order. It is convenient to define these 
expressions through the operators : X an (q) :. These are related to the rapidity operators as follows 

: R a , y (q) := R° an (q+ : X a „(q) :) - < 7 (q) , (96) 

where i?i° 7 (g) is the suitable GGS eigenvalue of R a ^{q) [see Eq. (73)]. 

The operators : X a ^(q) : contain the same information as the rapidity operators, and involve, exclusively, the two- 
pseudoparticle phase shifts ^ a ,-y-a' ,^'(<]j l') we define below. In Appendix B we introduce (96) in the BA equations 
(66) — (68) and expand in the scattering order. This leads to 

oo 

:i a)7 (?):=^iW(9), (97) 
j'=i 

where j gives the scattering order of the operator term Xa)y(q)- For the first-order term we find 

q' ,a' ,7' 

where the phase-shift expressions are given below. We emphasize that while the expressions for the phase shifts 
$ ctj7 . Q ,/ i7 /((7, q') are specific to each model because they involve the spectral parameters, the form of the operator term 

Xa)y{q) (98) is universal and refers to all the solvable electronic multicomponent quantum liquids. 
It is useful to introduce the phase shifts ^ a ^- a ',-y' such that 

~ , - / sinfc sinfc'\ , . 

$c,0;c,o(fc, k ) = $ c ,0;c,0 , , (99) 

\ u u J 



' sin k 



(100) 



*Q, 7 ;c,o(r, k') = $ Q , 7 :c,0 



sin k 1 ~ 



(101) 



*a, 7 ;a, 7 ' {r, T J ) = $a,-y;a, 7 ' (r, r') , (102) 

where 7 > and the phase shifts $ a , 7;Q ' !7 ' are defined by the integral equations (B30)-(B40) of Appendix B. 
The two-pseudparticle phase shifts can be defined in terms of the phase shifts $ Q . 7;a ', 7 ' as follows 

$a, 7 ;a',y (<?,</) = $a, 7 ;a', 7 < (^(ff).^^) • (103) 

The quantity ^> a , 7 ;Q' , 7 '(<?j l') represents the shift in the phase of the a', 7' pseudoparticle of pseudomomcntum q' due 
to a zero- momentum forward-scattering collision with the a, 7 pseudoparticle of pseudomomentum q. 

In Appendix C we use the Hamiltonian expression (83) in terms of the rapidity operators to derive the expression for 
the normal-ordered Hamiltonian. We find that in normal order relative to the suitable GGS (or GS) the Hamiltonian 
Hso{i)i Eq- (83), can be written as 
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■■H so{4) :=J2H {3) ■ (104) 

For example, the first and second pseudoparticle-scattering order terms read 

d (1) =E^) : ^):. (105) 

and 

^ (2) = ^ E E 7|/«,7;«',7'(<?,9') : ^a,7(3) ■■ N a ,„,(q') :, (106) 

a 9,a,7 q',a',j' 

respectively. All the remaining higher-order operator terms of expressions (97) and (104) can be obtained by combining 
the rapidity equations (66) — (68) and the Hamiltonian expression (83). 

In Appendix C it is shown that the pseudoparticle bands e" (g) can either be expressed in terms of the phase shifts 
(103), with the result 



°o(«) = ---2tcosK(°\q) + 2t J ^ ^ dk<5> cfi , cfl (k, (g)J sin k , (107) 



; „Q<+) 

< 7 (<?) = -7^ + 4iiieVl-« 2 [4^(9)-H] 2 + 2i y ( ) dfc$ c , 0;c , 7 (fc, sinfc 



for 7 > and 

'q(-) 

for 7 = 0, 1, 2, or simply as 



(108) 



e° 7 (g) = 2< y ( ; dfe*c,0;», 7 (fc,i2^ 7 (g)J sinfc, (109) 



c°o(«) = e'oW + / dk2tr 1c , (k) , (110) 

e° 7 (g) = / dr2t Vcn (r), (111) 



for 7 > and 

/■<>(<?) 

e° 7 (g) = / dr2^ S;7 (r), (112) 

*/ 00 

for 7 = 0, 1, 2, where the fuctions 2tr] C fl(k) and 2t-q an (r) are the solutions of the integral equations (C7)-(C9) of 
Appendix C. The GGS (or GS) rapidity functions K^\q) and i?i° 7 (g) are defined by Eqs. (77), (78) combined with 
Eqs. (B4)-(B6) of Appendix B. Combining Eqs. (70) and (71) with the band expressions we find that the bands 
(108), (111) and (109), (112) vanish at the limits of the pseudo-Brillouin zones, i.e. 

4, 7 (9$) = 0. (113) 
The associate group velocities, v an (q), are given by 

Va «ti) = — dq ' ^ ^ 

The two-pseudoparticle forward-scattering phase shifts $ Q , 7;Q ', 7 ' (g, q') defined by Eq. (103) and velocities 
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«q, 7 = v an (q Fa ,y) , (115) 

play an important role in the physical quantities when criticality is appoached |p5| , p6f . 

All Xa)y(q) terms of the rhs of Eq. (97) are such that both the / functions of the rhs of Eq. (106) and all the 
remaining higher order coefficients associated with the operators H^' of order j > 1 have universal forms in terms of 
the two-pseudoparticle phase shifts and pseudomomentum derivatives of the bands and coefficients of order < j. This 
follows from the fact that the 5*-matrix for j-pseudoparticle scattering factorizes into two-pseudoparticle scattering 
matrices, as in the case of the usual BA S- matrix ||, 3^|?9|,[l0|] . For example, we show in Appendix C that although 
the second-order term of Eq. (106) involves an integral over the second-order function Xa)r(q) [ see Eq. (C14) 
of Appendix C], this function is such that can be written exclusively in terms of the first-order functions (98), 

q,a.~f a, 7 j — ±1 

where |v a 7 | = C a . ) = 1 for x > 0, and 9(x) — for x < and the "Landau" / functions, f an - 7 ', a '(q,q'), 

are found in that Appendix to have universal form which reads 

fa 27rt7 Q;7 (g)$ Q;7;a ' 57 '(g , ; q') + 2nVa> n > (q')<& a f^/. a ^(q f , q) 

oo 

+ 2?r E E E ^(^a'^7'0^c t '',7'' ^; c t '',7''*a'^7'';a,7(J9Fa'^7'',g)^ , Q '^7''; Q ^7'(jgFa'^7'',g') , (117) 
j — ±1 a" 7"— 

where the pseudoparticle group velocities are given by Eqs. (114) and (115). 

We can also write the expression of for higher scattering orders j > 2. The main feature is that for all sub- 
canonical ensembles, energy scales, and pseudoparticle scattering orders, the Landau-liquid terms of that Hamiltonian 
involve only zero-momentum forward-scattering. 

Let us consider that the normal-ordering of the Hamiltonian (104) refers to the initial GS. In this case, and in 
contrast to the low-energy Landau theory studied in Refs. p0| , p2] , ^3| , p4||3^ ] , the Hamiltonian (104) describes finite- 
energy transitions involving a small density of excited pseudoparticles. Furthermore, a finite-size analysis |26| confirms 
that the excitations which control the response and transport properties at finite energies involve a small density of 
excited pseudoparticles. Fortunately, their description is thus within the range of our generalized pseudoparticle 
perturbation theory. It is convenient to define the normal ordering relative to the initial GS but with energy shifted 
by 

uj Q = lim [E GGS - E G s] • (118) 

This leads to a small (to — loq) energy description. Fortunately, the relevant GS - GGS transitions correspond to 
energy values ujq associated with small-pseudoparticle-density excitations J2f||2(|. For each of these loq values there 
is a low-energy (u> — u>o) pseudoparticle-perturbation theory. The associate Hilbert subspace is spanned by the set 
of GGS's of energy w and by states originated from these by pseudoparticle-pseudohole processes. The latter states 
have small momentum and low (lo — loq) energy relatively to the initial GS (we have shifted its energy by loq). Normal- 
ordering relative to the energy-wo shifted initial GS is equivalent to separate the problem (104) into (a) the GS - GGS 
transition of finite-energy (118) and (b) a low (u) — loq) energy theory associated with both the corrections of order 
1/N a to the GS - GGS-transition energy (118) and the Landau-liquid pseudoparticle-pseudohole processes around the 
final GGS. In most cases the excitation energy (118) associated with the GS - GGS transitions is finite. The finite 
energies are the energy gaps of the states II and (or) non-LWS's and non-HWS's relative to the initial GS. 

At finite values of S° and the low-energy Hilbert space is entirely spanned by states I |9| p0|j36] | . We note that 
for these states fixing the electron numbers fixes the pseudohole numbers fl20fl and the a, 7 bands are empty for 7 > 0. 
Therefore, if we fix the electron numbers we have that at energy scales smaller than the gaps (118) for the states 
II and non-LWS's and non-HWS's only Landau-liquid excitations (b) are allowed. This justifies the Landau-liquid 
character of the problem at low energies |Q,||j|j2rj. Obviously, if we change the electron numbers, there occur at low 
energy GS - GS transitions which are particular cases of the above GS - GGS transitions (a). Understanding these 
GS - GS transitions permits one to express the electron in terms of pseudoholes and topological momentum shifts 

Let us evaluate the general luq expression for GS - GGS transitions such that the final sub-canonical ensemble is 
characterized by vanishing or small values of N h _ la /N^. Here l a = ±1 is defined by Eq. (7). We also assume that 

a,— — 

both the initial GS and final GGS correspond to canonical ensembles such that 5" 7^ and belonging the same sector 
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of parameter space, i.e. with the same l a numbers. We note that since GS's are always S a LWS's I or HWS's I there 
are neither a, ^If- pseudoholcs [sec Eqs. (37) and (38)] nor a, 7 heavy pseudoparticles (7 > 0) in the initial GS and 
that N* = [see Eq. (13)] for that state. From Eqs. (104) - (106), we find for the GS - GGS gap 

^0 = E e%(q Fafi )AN afl + £„(0)N, n + 2fiAS c z + 2fi HAS s 2 , (119) 

a 7=1 

where the band expressions at the pseudo-Fermi points are given in Eq. (C12) of Appendix C. Inserting in Eqs. (11) 
both the pseudohole expressions (12) and the expressions S z = — h[N a — N] and S z = — i [iV^ — TV J we find after use 
of Eqs. (13) and (20) the following general expressions for AN C £ and AN S $ 

AN cfi = -2l c AS c 7 - 2N* - V 2 7 iV Ci7 = - [N h _, c - AN h lc } - 2N Z C - V 2 7 iV Ci7 , (120) 

7=1 7=1 

and 

a 7=1 7=1 

= ~\ Et<^ " A7V a% + ^ " E ^c )7 ^ E(! + i) N >« ■ ( 121 ) 
a 7=1 7=1 

Generalization of the chemical-potential and magnetic-field expressions (47) of Ref. |33| to all sectores of parameter 
space leads to 

H = -<o(gFc,o) ~ es > o(q 2 Fsfi) , 2 Mo | J ff|=-e 5 %(g Fs , ). (122) 

From the use of Eqs. (20) and (120) — (122) in Eq. (119) leads to the following gap expression 

oo 

= 2\fi\N^ +2^\H\[N^_^ + £jV, l7 ] +E<7(°)^.7> (123) 

7=1 7=1 

which is the general expression for the gap for GS - GGS transitions. It follows from Eqs. (13), (19), and (20) that 
the gap expression (123) can be rewritten as 

u, = 2| M | (E 7^,7 + K) +2^\H\ (£(1 + 7)A^7 + K) + ]T <? a „{W;i • ( 124 ) 

7=1 7=1 7=1 

In reference |2(| the concept of GGS is generalized to a, 7 > pseudoparticle filled seas with pseudo-Fermi points 
either of the form (a) qFa.-y — T ^°'] v Na -^ r (b) qFa.-y = 7rj ^ a ' 7 ■ We note that our definition (43) corresponds to the 
case (a) for a = c and (b) for a = s. While the first two terms of the rhs of Eq. (124) are the same for either case, 
the third term vanishes for the case (a) [this follows from Eq. (113)] and reads e° (0) for a, 7 seas corresponding to 
the case (b). This explains the form of the latter term in Eq. (124) which refers only to the s,7 pseudoparticles. 

For GS - GGS transitions to pure LWS's or HWS's II the use of Eq. (13) leads to 

N* = 0, (125) 

for both a — c,s. It follows that for such transitions the gap is given by expression (124) with N* = and N* = 0. 
Therefore, the single (7 > 0) a, 7-heavy-pseudoparticle gap is 2-f\/i\ and 2(l+7)/io|i?| +e" 7 (0) for c and s, respectively. 

The gap expression (123) — (124) involves the chemical potential magnetic field H, and the zero-momentum a, 7 
(with 7 > 0) band values. For constant values of U the above quantities change as follows: For fixed values of S z , 
\fi\ changes between = A " H for — * (and n — ► 1) and a maximum value, \fi\ = Ar ^ H + &t, for S° — > ±^ 
(and both n — ► and n — > 2), where Amh is the half-filled Mott- Hubbard gap jl8],[52]]. The gap Amh is both an 
increasing function of U and of \S S Z \ [[32|. On the other hand, for fixed values of 5f! 7jH\ changes between \H\ = 

for jf > (and m — > 0) and a maximum value, \H\ = H Cl for S s z —> ±-y (and m — > ±n and n — * 2), where iJ c is 

the critical field for onset of fully polarized ferromagnetism defined in Refs. (3^j34|]. [Its expression presented in these 
papers refers to the parameter-space sectors where the electronic density is such that < n < 1. The generalization 
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of that expression to the sectors where the density is such that 1 < n < 2 is obtained by replacing in it n by 2 — n.] 
Finally, for fixed values of S c z the band e° s 7 (0) changes between e° (0) = for jf- — ► and a minimum negative value. 
A more detailed study of the U, density, and magnetization dependence of the gap expression (123) — (124) will be 
presented elsewhere. 

Given an initial GS and a final sub-canonical ensemble we can write 

: H \= H + HLandau , (126) 

where Hq has eigenvalue loq and corresponds to the GS - GGS energy (a) and HLandau is normal-ordered relative to 
the initial GS (with energy shifted by u)q) and is of the form 

HLandau = + , (127) 

where 

H { l = E e ^) : &vr(«) : . (128) 



ea,o(«) = 4,o(9) - 4>Fa,o) , e a „(q) = C^(q) - 6 a , s e^(0) , (129) 

* (2) 

and H L ' is given in Eq. (106). The Hamiltonian (127) describes both the GS - GGS transition and Landau-liquid 
pseudoparticle - pseudohole excitations relative to the GGS. (The c,0 and s,0 pseudoparticle bands are shown in 
Figs. 7 and 8, respectively, of Ref fl32fl .) Here (127) are the Landau-liquid Hamiltonian terms which are relevant at 
low energy (lo — ujq). Therefore, the corresponding second-order Hamiltonian is suitable to study the physics at low 
positive energy above the gap, i.e., small (li — ojq) p5| , p6| . 

In general, the different final GGS's of the Hilbcrt subspace where a given initial GS is transformed upon excitations 
involving a small density of pseudoparticles have different energies (123) — (124). This implies that the study of the 
quantum-liquid physics at energy-scale co involves, in general, transitions to one sub-cannonical ensemble only. When 
two or several possible final GGS's have the same energy gap (loq) the physics involves the Hilbert subspace spanned 
by all Hamiltonian eigenstates associated with the corresponding different final sub-canonical ensembles. 

For sub-canonical ensembles such that S" ^ and characterized by small values of /N„ [where l a is given 

in Eq. (7)], the gap expression (123) leads directly to 

oo 

Ho = 2HA^_ ¥ + 2^ \H\[N h s _ lf + + E e $(0)#a, 7 ■ (130) 

7—1 7—1 

There is a remarkable similarity between the general Hamiltonian (126) and its version in the Hilbert subspace 
spanned by the states I dj2rj,Q. Moreover, in the case that the GGS is the GS it self or a GS differing from it by 
one or a few electrons the gap (123) — (124) vanishes and the Hamiltonian (126) reduces to the above low-energy 
Hamiltonian. In this case it refers to gapless GS - GS transitions which change the electron and a,0 pseudoparticle 
numbers and to Landau-liquid a,0 pseudoparticle-pseudoholc excitations around the final GS. 

The two-pseudoparticle forward-scattering phase shifts ^ a ,T,a',f'(l> l') defined by Eqs. (99) — (103) and (B30)-(B40) 
and the velocities (115) play an important role for the physics at energy scales around ujo pq , |26| . [In the case of 
7 = 0, the velocities (115) are plotted in Fig. 9 of Ref. |32[ -1 

In this section we have presented the expression of the normal-ordered Hamiltonian and rapidity operators in the 
pseudoparticle basis. This confirms the consistency of the one-dimensional Landau-liquid theory which was shown 
to refer to this basis. The advantage of using the pseudoparticle basis is that (i) the problem becomes perturbative, 
i.e. from the point of view of the pseudoparticle interactions it is possible to classify which scatterings are relevant 
and (ii) we can describe the relevant high-energy GS transitions because these refer to a small density of excited 
pseudoparticles and are, therefore, within the range of the perturbation theory. 

For energies ui just above the set of energy values uj of the form given in Eqs. (123) — (124) only the two- 
pseudoparticle interactions are relevant, as can be confirmed by finite-size studies at the corresponding critical point 
|26| . Therefore, and for simplicity, we have presented in this paper only the two first Hamiltonian terms of expression 
(104). However, Eqs. (66) — (68) and (83) contain full information about all the remaining terms of higher-scattering 
order. All these many-pseudoparticle terms only describe zero-momentum forward-scattering interactions. 
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VI. CONCLUDING REMARKS 



In this paper we have introduced a pseudohole and heavy-pseudoparticle representation for all the Hamiltonian 
eigenstates of the Hubbard chain in a magnetic field and chemical potential which is valid for all sectors of Hamiltonian 
symmetry. In this picture all Hamiltonian eigenstates can be generated from a single reference vacuum, the half-filling 
and zero-magnetic-field GS. 

Our algebraic approach permits an operator analysis of the BA solution and goes beyond that solution by generating 
the non-LWS's and non-HWS's of the ?7-spin and spin algebras. We find that all GS excitations can be divided into (a) 
a GS - GGS topological transition which changes the a, (3 pseudohole and a, 7 heavy-pseudoparticle numbers and (b) 
one pseudoparticle-pseudohole excitation around the final GGS. It will be shown elsewhere |2^,^(| that the quantum- 
liquid physically important transitions are the GS - GGS transitions which occur at energies ojq of form (124). Further, 
for low (u> — wo) energies our algebraic approach leads to the construction of a GS normal-ordered Hamiltonian which 
has in the pseudoparticle basis a universal form involving only k — 0, forward-scattering pseudoparticle interaction 
terms. 

A central result of our paper is that the perturbative character of the pseudoparticle basis also refers to high- 
energy states provided that the density of excited a, 7 pseudoparticles is small. This perturbative character implies 
that the two-pseudoparticle Landau / functions and forward-scattering amplitudes are finite, in contrast to the non- 
perturbative electronic representation, in which the two-electron forward-scattering amplitudes and vertices diverge. 

The pseudoparticle algebra allows us to see immediately the origin of the "universal" character of the one- 
dimensional integrable quantum liquids at all energy scales, which can be understood as a straightforward gener- 
alization to pseudoparticles of Wilson's low-energy renormalization group arguments: the pseudo-Fermi points of the 
pseudoparticle GS's and GGS's replace the particle Fermi points, and close to the pseudo- Fermi points only few types 
of two-pseudoparticle scattering processes are relevant for the low (lj — u>o) energy physics. 

Concerning applications of our ideas to other theoretical models and to problems in real materials, we note that it 
is of considerable interest to determine the extent to which the algebraic structure and pseudoparticle perturbation 
theory remain valid in systems which are not integrable but which behave as electronic Luttinger liquids; an example of 
such a system is the one-dimensional extended Hubbard model. In terms of real materials, it is known that the present 
one-dimensional quantum systems provide useful (albeit idealized) descriptions of the physics of quasi-one-dimensional 
solids. 

Finally, we turn to directions for further research. In addition to some specific calculations mentioned above, our 
immediate aims are to apply the concepts and techniques developed here to study the finite-frequency spectral and 
transport properties of the Hubbard chain. The occurrence of only zero-momentum pseudoparticle forward scattering 
combined with the infinite pseudohole- and pseudoparticle-number conservations laws can be used to extract important 
information on the finite-frequency behavior of correlation functions ^6|. We will show elsewhere |25| that the a, 7 
pseudoparticles are the transport carriers at all energy scales and couple to external potentials and that the topological 
GS - GGS transitions introduced in this paper give rise to the finite-frequency-conductivity absorption edges. We 
emphasize that our pseudoparticle operator basis provides a more suitable starting point for future investigations on 
that open problem. This is a subject of considerable interest for the understanding of the mechanisms behind the 
unusual spectral and transport properties of the novel low-dimensional materials. 
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APPENDIX A: RELATION TO THE USUAL BA NOTATION 



Our pseudoparticle operator representation required a suitable choice of notation which is directly related to the 
BA notation used by Takahashi |19| . In this Appendix we relate the pseudoparticle notation to his notation. In our 
c, 7 and s, 7 parameters below we consider that 7 > and 7 = 0, 1, 2..., respectively. 

The N an pseudoparticle numbers are related to the numbers N, M', M' n , and M n of Ref. [^9| as 

JV Ci0 = N - 2M' , jV c , 7 = M' 7 , N.„ = M 1+7 . (Al) 

Moreover, our onsite parameter u = % equals the interaction U of that reference. Considering the discrete pseudo- 
momentum qj of Eq. (28) , our rapidity eigenvalues are related to the rapidity-real parts kj , A'™ , and A™ of Ref. Jl^] 
as follows 

A'] Aj +7 
K( qj ) = kj, R cn (q,) = -^, R s ^ qj )=^—_ (A2) 

Finally, the relation of the pseudoparticle quantum numbers of Eq. (28) to the Takahashi quantum numbers Ij, J' a , 
and J™ is 

r/ I,- /;•" ir^j) + ~<- (A3) 
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APPENDIX B: NORMAL-ORDERED RAPIDITY OPERATOR EXPRESSIONS 



Following the discussion of Sec. V, the perturbative character of the system implies the equivalence between 
expanding in the pscudoparticlc scattering order and/or in the pseudomomcntum deviations (95). In this Appendix 
we give a short description of the calculation of the normal-ordered operator expansion for the pscudoparticlc rapidities 
(97). We focus our study on the evaluation of the functions Xa^q) and Xa)y{q) and associate functions Qa)y(q) and 
Qa)y{q) (eigenvalues of the operators Xa \q) and Xa \q) and Qa\q) and Qa\q), respectively). 

We evaluate here the first-order and second-order terms of the eigenvalues of the operator (97). Equation (86) then 
allows the straightforward calculation of the corresponding operator expressions. 

In the thermodynamic limit, Eqs. (66) — (68) lead to the following equations 

m =5--£ /ff HN S M) tan- (^M_^lM) 



7=1 J 1> 



,(+) 



^djN^Wtn- ^M-^W ), (Bl) 



2Resm- 1 (u(R c .,(q)-i 1 )) =q+- f" dq'N cfi {q>) tan" 1 ^Wh 

V / 7T Jgi-) V 7 1 

i r q(+} ' 

+ T- E /,7 ^'^C7'('Z')e7,7'(^7('Z) - Rc^(q')) , (B2) 



and 



0=5-- , 

7T ./„(-) 



1 / a W )'~-' C- ( '': ft ' l ' ,) ) 

J g (J V 1 + 7 7 

i r qi+) > 

+ ^rJ2 /,7 ^'^,7'(9')0i+7,i+7'K,7(9) - R,,ytfj) ■ (B3) 

We start by considering the GGS eigenstate rapidities of Eqs. (72) and (73). If we insert in Eqs. (66) — (68) the 
GGS distribution (75) after some algebra we find that the functions of the rhs of Eqs. (77) — (78) defining the inverse 
of these GGS rapidities are solutions of the following integral equations 

cosfc f r<s -<> , 2irp s0 (r) . 
27rp c (k) = 1 + / dr— ? 1 ■ v ; , = , B4) 



27rp c>7 (r) = 2Re(—= 



i/l — ?i 2 [r — z 7 ] 2 ' JQ(-) 



7T7 



l + ( = ) 2 



(B5) 



and 



27rp s 



,(r) = / dfc- 



27rp Cj0 (A:) 



Q<-> 7r (i+ 7 )^i + (-^zi)2 

^ eW 1+7 (r-r') 

, , dr'— ^ / -2 7 r Ps ,o(r'), 

(-) Z7T 

s,0 



(B6) 



where 
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1 7-1 2 

= ^ 7 ' { 7[i + (^) 2 ] + § ^[i + (i) 2 ] } 



| 7 -7'|[l + (^) 2 ] (7 + 7')[l + (^y) 2 



2+2^__[2j_2iL _ i 



X — ^ ^ 

+ C7^1TWTT=^=H } ' (B7) 

is the derivative of the function (69) and the parameters and r^fy are defined combining Eqs. (77) — (78) and 
(79). We have omitted in the rhs of Eqs. (B4)-(B6) terms associated with contributions of order j > 2 in the density 
of excited pseudoparticles. These terms do not contribute to the quantities to be evaluated in this paper, as we discuss 
below. 

Let us now consider small deviations from a GGS or GS. The eigenvalue form of Eq. (96) is 

SK(q) = K^(q + 6X c ,o(q)) - K^(q) , (B8) 

for a = c and 7 = and 

5R an {q) = R° an (q + 8X a „(q)) - R^q) , (B9) 

for all remaining values of the quantum numbers a and 7. Here 8K(q), 8R a ^(q), and 8X a ^(q) are the eigenvalues of 
the operators : K{q) :, : R aa {q) :, and : X an (q) :, respectively. From Eq. (97) 8X an {q) can be written as 

5X an (q) = Xi)\(q) + X^ (q) + ... , (BIO) 
where Xa\(q) is the eigenvalue of the operator Xa\(q). Expanding the 8R an (q) expressions (B8) and (B9) we find 

tf(g) = f>«(«), (Bll) 

i=0 

and 

00 

R a , 7 (q) = J2^( ( l)' (B12) 

i=0 

respectively, [and 8K{q) = J^iti^^Hl) an d 8R a!j (q) = ^0,7(9)] where the zero-order GGS functions JsT^^g) 

and i?i°^(g) arc defined by Eqs. (77) — (78). From the resulting equations we can obtain all derivatives of the 

GGS functions K^(q) and Ra.ly(q) with respect to q. The terms of the rhs of Eqs. (Bll) and (B12) involve these 
derivatives. For instance, the first-order and second-order terms read 

K^q)= d J^X^q), (B13) 
R^ {q) ^ d -^Mx^{q), (B14) 



and 



K^ = d ^X^ qH l^ { xi% )? , (B15) 
R W {q) = ^M X m {q) + (,)] 2 , (B16) 
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respectively, and involve the first and second derivatives. 

From Eqs. (B1)-(B3) [with N an (q') given by the GGS distribution (75)] we find that the first derivatives dg f g >(<?) 

dit' ' (q) 

and — ^ — can be expressed in terms of the functions (B4)-(B6) as follows 



dK(°\q) _ 1 



dq " 2irp c , (KW{q)) ' 

and 

dR { °Uq) = 1 

dq " 2irp cn (R { °l ( (q)) ' 

The second derivatives d K ^ q l^ and d R ^ q ^ then read 



(B17) 



(B18) 



d 2 KW(g) 1 d Pe ,o(fc) 

<V - - [2*p c , Q {Km{qW 2 *^— k = K{a) ^ ' (B19) 

and 

By introducing both the distributions 

N an {q) = NM(q) + SN ar/ (q) , (B21) 

and the first-order and second-order functions (B13)-(B16) into Eqs. (B1)-(B3), we find after expanding to second 
order that for j = 1 and j = 2 the functions Xa)y (q) can be written as 

X^(q)=Q^(q)+Y^(q), (B22) 

where Ya,](q) i s even i R 1 

Y$(Q)=Y<${-q), (B23) 



and does not contribute to the physical quantities to second scattering order (and to second order in the density of 
excited pseudoparticles) . Up to this order only Qa\(q) contributes. 

To derive this result we have expanded the expression for the even functions Y"if 7 (q) for j = 1 and j = 2 to second 
order in the density of excited pseudoparticles. Following the pcrturbative character of the quantum liquid in the 
pseudoparticle basis, the obtained expression is exact up to j = 2 pseudoparticle scattering order. Moreover, we find 
that the even function vi j7 (q) does not contribute to the physical quantities up to that order and can, therefore, be 
omitted. 

Let us introduce the functions Qajy{r) such that 

QL 1 , ) 7 (5)=0i 1 ; ) 7 « ) 7 (?)), (B24) 
for all values of a and 7. Note that following Eq. (64) 

m (q) = ^M. (B25 ) 



u 



It is also useful to define the function Q^\k) such that 



Q^\k) = Q^A- (B26) 
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By introducing both the distributions (B21) and the first-order functions defined in Eq. (B24) into Eqs. (B1)-(B3), 
we find after expanding to first order that the functions Qa)y{r) are defined by the following system of coupled integral 
equations 



n J Q 



-E 



,„H 

,(+) 



is,-, 2 , / r — 

dqSN s ^(q) — tan 

7T VI 



,„(-) 



1 / 

dqSN Cyl (q)— tan 

7T V 



(0) 



+ 7 / 

,/r-i?^(g) 



7 



<-) 7r[l + (r-r') 2 



(B27) 



and 



r s w , 0( „i ta „-.(^M) 



+ 



7T7 



1 + ( r__2_)2 



«sw-rf'**^(,)i ta -(i^) 

^c,0 

/•9 (+> 1 /•'' <+) 

7 '=o- / 'i.y J <o 7t(1+ 7 ) [l + (^.)2 



r — r 



(B28) 



(B29) 



[The use of Eq. (B27) in Eqs. (B28) and (B29) allows the expression of both Qi]\(r) and Qi]\(r) in terms of 
free terms and integrals involving Oi^oV)-] Combining Eqs. (B27)-(B29) with Eqs. (B24) and (B26) leads to Eq. 
(98) with the phase shifts defined below. Note that at first order we can either consider the function Xa)y(q) or the 
associate function Qa)y(q) of Eq. (B22). Both functions are of the form (98). However, the general expressions for the 
phase shifts associated with the function Xa,\(q) of Eq. (98) have extra terms. These arise from the function Ya]](q) 
of the rhs of Eq. (B22). If we expand the physical quantities involving the phase shifts to first (and second) order in 
the density of excited pseudoparticles these extra terms lead to vanishing contributions. For instance, expression (98) 
is identical if we use in it either choice for the phase shift expressions. For simplicity, we omit here the phase-shift 
extra terms associated with the function Y^{q). We find that the phase shifts associated with the function Qajy(q) 
arc defined by the following coupled integral equations 



1 r r ^o $ n 
<&c,0;c,o(r,r') = - / dr"^ 

7T J r (- 1 + 



;co(rV) 



+ (r - r") 2 



(B30) 



1 r — r' 1 
* c ,o :c ,7 {r, r') = tan _1 ( ) + - 

7T 7 7T 



r iV $ n ( r " r') 
(-) T l + (r-r") 2 



(B31) 



1 , r — r' 1 f r ^° $ 
$c,o ;s> 7(r-,r / ) = tan-^— — ) + - / dr" - 

TT 1+7 7T J r (~) 1 



, // # a ,o;a, 7 (r",r') 



+ (r - r") 



ll\2 ' 



(B32) 



1 r — r' 1 ( Tc 
®c. T .c.o (r, r') = - tan"^ ) - - / dr" 



c,0;c,0 



(r",r') 



7 



7[l + (^) 2 ] 



(B33) 
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4>, 



(r,rO = ^(r-^)-ijf ( ; 



1 r r c° $ n / (r" r') 

I ^11 Y c,0:c,7' \i i i } 



7[l + (^-) 2 



(B34) 



*c, 7;s ,7' (r,r') =-- 



.(+) 



(B35) 



1 r — r' 1 

$s, 7 ;c,o (r, r') = tan _1 (— — ) + - 

7T ii- 



dr 



„ ®c,0;c;Q (r",r') 



1 + 7' nj r ij (i+ 7 )[i + (!^)2] 



(-) 

a,0 



dr"$ s>0]cfi (r",r') 



2tt 



(B36) 



1 r r c,o $ n , ( r " r >) 



W r H (l + 7 )[l + (^)2] 



e 



[i] 



r ( + ) 

s,o _ 1+7 

s,0 



i(r-r") 



2tt 



(B37) 



2tt 



„(+) 



0i+7,i+ 7 ' ( r - r 'j | /-^.o' $ n ,(r" r 1 ) 



„(+) 



e 



[i] 



^1+7 



2 it 



(B38) 



For 7 = Eqs. (B36) and (B38) can be rewritten as 



i r# 

$*o(^') = — tan" 1 (r-r') + / dr"G(r, r")$ s ,Q; C ,o (r", r') , 



(B39) 



and 



„(+) 



1 _i,r-r\ 1 /■ p ».» ,„ tan- 1 (r" - r') 

* a ,o 5 .,o r,r' --tan 1 — — - — / d r " — 

7r 2 7r^ 7 r (-) 1 + (r — r 'J^ 



,(+) 

,0 



+ / dr"G(r,r")$ s ,0;s,o(r"y) , 



(B40) 



and the kernel G(r,r') reads 



G(r,/) = -- 



1 + ((r - r')/2) 2 



1 z' . . , ,, l(r)-l(r') 
l-~lt(r)+t(r')+ y! r _ r \ 



(B41) 



with 



1 r 



<(r) = - tan _1 (r + r£t ) ) - tan _1 (r + r£ 

7T 



(B42) 



and 



l{r) = - [ln(l + (r + r«) 2 ) - ln(l + (r + r^) 2 ) 

TT I ' . 



(B43) 
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In order to evaluate the second-order functions Q^)y(q) of the rhs of Eq. (B22) for j = 2 [see also Eq. (BIO)] we 
introduce the functions (B13)-(B16) and distributions (B21) in Eqs. (B1)-(B3). Expanding to second order we find 
after some algebra 



Qi%(q)-Qi%*\q) + ~[[Q ( ^ 2 h 



where 



It is also useful to define the function Q^ 2 ^'{k) such that 

The functions Qa^\r) are defined by the following system of coupled integral equations 



~f=o J ^,-> ZnPs^yRsAiq) J ^(1 + 7 ) 



1 + ( z^gM )2 



s / : 



dqSN {q) _&M i 



+ E 



E 



[Q£W S , 7 )] 2 



o 27rp 5 , 7 (r 5 , 7 ) j=±1 2tt(1 +7) [l + C^ 1 ^) 2 



E 



W, 7 ) 



E 



[Q$(j<ZF C , 7 )] 5 



+ 



=1 2np cn [r cn j j=±1 2irj 
[1 + (r - r') 2 ] ' 



/ dr'- 

J r (-) 7T 



(B44) 
(B45) 

(B46) 



(B47) 



and 



Q&Hr) = -J*° dqSN cfi (q) 



cos K(°\q) Q$(<z) 



27rp c , (fc(°)(9)) 



7T7 



l + ( 



„(+) 



Y — 1 c.V 



7' 
COS 



E 



2irp c ,y(RW(qj) ^ 



(Q) j=±i27r 7 [l + (- 



3^,0 \2 
7 ' 



0(JV c , 7 ,) 1 



^EW?yOW)] a ^e^(r-jrcy) 



(2 '*V) 



7T7 



1 + (1=^)2 



„< + ) 



Ql,y(r)= /, , dqSN cfi (q) 

" , 9c,o " 



27rp c ,o(fc(°)(g)) ^(1+^ 



l + (^) 2 



(B48) 
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r 



+ 



COSQ 1 x ^ [QcfliilFcfl)] 2 



2^ c ,o(q) ,£±i 2n(l + 7) [l + C-^) 2 



+ 




s,0 



(B49) 



As for the first-order case, the use of either the function Q ( a,j(q) defined by Eqs. (B44)-(B49) or of the full function 

(2) 

X£!y(q) of Eq. (B22) for j = 2 leads to the same results for the physical quantities up to second order in the density 
of excited pseudoparticles. 

Note that the free terms of Eqs. (B47)-(B49) involve the first-order functions only. This implies that the unique 
solutions of these integral equations can be expressed in terms of the first-order functions. Therefore, following Eqs. 
(98) and (B22) the second-order functions can also be expressed in terms of the pseudoparticle phase shifts. 
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APPENDIX C: NORMAL-ORDERED HAMILTONIAN EXPRESSION 



In order to derive the first-order and second-order Hamiltonian terms of Eqs. (105) and (106), we again consider 
eigenvalues and deviations. The energy associated with the Hamiltonian (83) reads 

it n r q c + o 4tN r q( <^j i 

E SO( A) = -^f N cfl { q ) C osK(q) + ^Y. N c , 1 (q)Re^l-[u{R c , 1 {q)- ll )Y 

+ ^[f-% £ -E^]- (CD 

7 = 1 

To caculate the bands (107) — (109) we introduce in this energy Eqs. (B13), (BI4), and (B2I) and expand the obtained 
expression to first order in the deviations with the result 



A 4o(4) = T^if*' SN Ci0 {q)[-2tcoaK(-°\q) 

^c,0 
„( + ) 



s / : 

/>r K ' poo , 2r ■ 1 

/ dkQW(k)sink-4tJ29(Nc,,)[ ' + ]drQ^{r)Re( [ ~ 7j , )}, (C2) 



+ 5N Ctl {q)[AtRe^l- [u(R^>(q) - * 7 )] 2 - lU] 

7 

r Q {+) r A~J r°° / ,,2 



where the functions Q^\k) and Qi]\(r) are defined by Eqs (B26)-(B28). As in the case of the second term of the rhs 
of Eq. (B22), up to second scattering order (and to second order in the density of excited pseudoparticles) the last 
term of Eq. (C2) does not contribute to the energy expression. Therefore, it can be omitted. The use of Eqs. (B22), 
(B24), (B26), and (98) in (C2) then leads after some straightforward algebra to 

A 4o(4)= £< 7 (<7)<W a , 7 (9), (C3) 

with the bands given by Eqs. (107) - (109) in terms of the phase shifts (B30)-(B40). 

To obtain the equivalent band expressions (110) — (112) requires introducing in the energy expression (C2) the 
functions defined by Eqs. (B26)-(B29) and performing some integrations by using symmetry properties of the obtained 
integral-equation kernels. After some algebra we find again the energy expression (C3) with the pseudoparticlc bands 
given by 



c 



\(q) =~- 2tcosK(°\q) - jf*"° dr^oW^tan- 1 (r - R^(q)) , (C4) 



e° -yO?) = -lU + UReJl - u^R^iq) - ry] 2 - / dk- tan 
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e° 7 (g) = -2t / 

JQ 



Q <+) ! 

dk— tan 1 

Q<-) n 



„<+) 
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1 + 7 
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<+) tan 



2tt 
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J r (-) 7T 



1+7 



2 [l + (r-r') 2 ] 



(C5) 



(C6) 



where the functions 2trj c ^ (k) and 2ti] a ^(r) are defined by the integral equations 
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2tr]c,o(k) = 2isinfc + 



cos A: 



,r (+) 

/ dr 



2*»7a,o(r) 



7T [1 + (r - ^) 2 



(C7) 



2tr/ C:7 (r) = -4tRe( 



yjl — u 2 [r — i^} 2 



2ir/ Cj0 (/c) 
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Q(") tt 7 



, eW +7 (r'-r) 
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s,0 



(C8) 



2tt 



(C9) 



The use of these equations and comparision of the band expressions (C4)-(C6) with the derivatives of the functions 
(C7)-(C9) readly leads to the simple expressions (110) — (112). 

In order to derive the expression for the second-order Hamiltonian (106) and associate / functions (117) we expand 
the energy (CI) to second-order with the result 



AE 



N a 
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,( + ) 
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dq5N cfi (q) , ... \ Qc,0 («) 



27rp c ,o(^(°)(9)) 



+ 



E 

=i • 

2isinQ 1 
2irp cfi (Q)2 



2rn(0)/ \ ■ I 



-AtRei 



dq5N Cyl (q)— , , . 



E muiFco) 



j=±i 



2 it±i 2 ^c, 7 (r c , 7 ) Vi_ u 2 [rci7 _i 7 ]2 



+ E 

7=1 

+ / dkQ {2 '*\k)2tsmk 

r(-) 



r) 



+ / wg( 2 7 *)(r)W— 

-oo "M^ x V 1 



m 2 [r — 17] 



■ u 2 [r — i"f]' 
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(CIO) 



Again, the last term of the rhs of this equation does not contribute to the energy expression to second scattering 
order (and to second order in the density of excited pscudoparticlcs) . Inserting the suitable functions in the rhs of Eq. 
(C14), performing some integrations by using symmetry properties of the kernels of the integral equation (B27)-(B29) 
and (B47)-(B49), and replacing deviations by pseudomomentum normal-ordered operators (84) we find after some 
algebra expression (116). Note that replacing in Eq. (116) the function Xa)y{q) by the associate function Qa}j(q) 
leads to the same result. By the use of Eq. (98) expression (116) can be rewritten in terms of the / functions (117) 
as given in the rhs of Eq. (106). 
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